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Abstract 

We review in this article the influence of surface waves on the thermally excited 
electromagnetic field. We study in particular the field emitted at subwalength dis- 
tances of material surfaces. After reviewing the main properties of surface waves, we 
introduce the fluctuation-dissipation theorem that allows to model the fluctuating 
electromagnetic fields. We then analyse the contribution of these waves in a variety 
of phenomena. They give a leading contribution to the density of electromagnetic 
states, they produce both temporal coherence and spatial coherence in the near 
field of planar thermal sources. They can be used to modify radiative properties 
of surfaces and to design partially spatially coherent sources. Finally, we discuss 
the role of surface waves in the radiative heat transfer and the theory of dispersion 
forces at the subwavelength scale. 
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1 Introduction 



Many condensed matter properties are determined by surface properties. Very 
often, surface waves, which are electromagnetic eigenmodes of the surface, 
play a key role. Let us mention a few examples. It has been demonstrated [1] 
that the lifetime of a molecule varies dramatically when a metallic surface is 
brought at a distance smaller than a micron. This effect is due to the resonant 
excitation of surface plasmons. The van der Waals force between a molecule 
and an interface is proportional to l/|e + Ip where e is the dielectric constant 
of the medium. There is therefore a resonance for the particular frequency 
such that e = — 1. This condition coincides with a branch of the dispersion 
relation of a surface wave. It can be viewed as a resonant excitation of surface 
charge oscillations. It was shown in [2] that the van der Waals force between a 
molecule and a surface can become repulsive depending on the relative position 
of the molecule and the surface resonances. Enhanced scattering due to the 
resonant excitation of surface charges has also been demonstrated for SiC in 
the infrared : a tip brought close to a surface generates a very strong scattering 
signal for some particular frequencies corresponding to the excitation of surface 
waves [3] . Both experiments can be understood by replacing the interface by an 
image whose amplitude is very large owing to the excitation of a resonance of 
the surface charges. Surface Enhanced Raman Scattering (SERS) is partially 
due to the enhancement of the electromagnetic field at the interface due to 
the excitation of a surface wave. The resonance of the electromagnetic (EM) 
field at an interface is also responsible for the enhanced transmission of a 
metallic film with a periodic array of holes [4,5]. The resonance of the EM 
field associated with the surface mode is responsible for the so-called "perfect 
lens" effect [6]. A key feature of all the above examples is that they involve 
the interaction of a surface and an object in the near field of the structure. As 
it will be explained in details in Section 2, surface waves are evanescent waves 
whose amplitude decreases away from the interface on a wavelength scale. In 
the far field, the influence of such modes is therefore negligible. In the near 
field on the contrary, their role is essential. 

We will see in Section 3 that surface waves can be excited by thermal fluc- 
tuations inside a body. The role of surface waves in the modification of the 
density of EM states at the interface has a strong influence on the thermally 
emitted flelds. Their intensity is many orders of magnitude larger in the near 
fleld than in the far fleld [7]. In addition, they are quasi monochromatic in 
the vicinity of the surface. This entails that their coherence properties are ex- 
tremely different from those of the blackbody radiation [8]. There have been 
recently several experiments that have probed these thermal flelds in the near- 
fleld regime : heating of trapped atoms [9], realization of a spatially partially 
coherent thermal source [10]. After reviewing these experiments, we will show 
how an EM approach with random fluctuating thermal sources can be used 
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to describe and analyse these effects. It is based on the fiuctuation-dissipation 
theorem. Wc will see that the knowledge of the electromagnetic energy den- 
sity gives acces to a fundamental concept : the local density of EM states. In 
section 4, we study the EM coherence properties near a material supporting 
surface waves and held at a temperature T. We will see that the emitted field 
has very peculiar spatial coherent properties in the near field. Indeed, the field 
can be spatially coherent over a length larger than several tens of wavelength. 
We then use this property to design coherent thermal sources. In section 5 
and 6 we show that the radiative heat transfer is enhanced by several orders 
in magnitude in the near field when two material supporting surface waves are 
put face to face. We will consider three cases : two nanoparticles face to face, a 
nanoparticle near a plane interface and two semi-infinite half-spaces separated 
by a narrow gap. In the last section, we will analyse the role played by the 
surface waves in the Casimir force i.e in the force of interaction between two 
semi-infinite bodies. We will see that this force is dominated in the near field 
by the interaction between surface waves. Finally, we review the work done 
to analyse the contribution of fiuctuating electromagnetic fields to the friction 
forces. 



2 Introduction to surface electromagnetic waves 



In this section, we give a brief introduction to the main properties of elec- 
tromagnetic surface waves. This particular type of waves exists at the inter- 
face between two different media. An electromagnetic surface wave propagates 
along the interface and decreases exponentially in the perpendicular direction. 
Surface waves due to a coupling between the electromagnetic field and a reso- 
nant polarization oscillation in the material are called surface polaritons. Prom 
a microscopic point of view, the surface waves at the interface of a metal is a 
charge density wave or plasmon. It is therefore called surface-plasmon polari- 
ton. At the interface of a dielectric, the surface wave is due to the coupling 
of an optical phonon with the electromagnetic field. It is thus called surface- 
phonon polariton. Plasmon polaritons and phonon polaritons can also exist 
in the whole volume of the material and are called polaritons. More details 
about this subject can be found in textbooks such as Kittel [11], Ashcroft and 
Mermin [12] and Ziman [13]. In what follows, we will focus our attention on 
surface polaritons propagating along a plane interface. Excellent reviews of 
the subject can be found in [14,15,16,17]. 
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2.1 Surface polaritons 



Let us now study the existence and the behaviour of surface polaritons in 
the case of a plane interface separating two linear, homogeneous and isotropic 
media with different dielectric constants. The system considered is depicted 
in Fig.l. 




Fig. 1. A plane interface separating medium 1 (dielectric constant ei, magnetic 
constant /xi) and medium 2 (dielectric constant 62, magnetic constant ^2) 



The medium 1 (dielectric constant ei and magnetic constant /xi) fills the up- 
per half-space 2; > whereas medium 2 (dielectric constant 62 and magnetic 
constant 112) fills the lower half-space 2; < 0. The two media are supposed to be 
local and dispersive so that their complex dielectric and magnetic constants 
only depend on cu. 

The three directions x.y.z shown in Fig.l arc characterized by their unit 
vectors x, y , z. A point in space will be denoted r = {x, y, z) = ,xx + yy + zz ~ 
(R, z) where R = -|- yy. Similarly, a wave vector k = {k^, ky, k^) will be 
denoted by k = (K, 7) where K is the component parallel to the interface and 
'y — kz the component in the z direction. 

A surface wave is a particular solution of Maxwell's equations which prop- 
agates along the interface and decreases exponentially in the perpendicular 
directions. Because of the translational invariance of the system, it can be 
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cast in the form: 



Ei(r,c^) 



^ Ex,2 ^ 



exp[i(K.R + 7i^)] (Mediuml), 



(1) 



E2(r,u;) = 



Ey,2 



exp[i(K.R - 722;)] (Medium2), 



(2) 



where 71 and 72 are given by, 



ei/ii/cQ — with /m(7i) > 0, 

72 = €2/^2^:0 - K"^ with /m(72) > 0. 



(3) 
(4) 



Here — uj/c where c is the speed of hght in vacuum. We now look for the 
existence of surface waves in s (TE) or p (TM) polarization. In what follows, 
we shall assume that the wave propagates along the y-axis. 



2.1.1 s -polarization (TE) 

In s-polarisation, the electric field is perpendicular to the plane {y,z). The 
electric field E is thus parallel to the x direction 



Ei(r,u;) = £;^,ixexp[i(K.R + 7i^)], 
E2(r,a;) = £;^,2xexp[i(K.R - 72^)]. 



(5) 
(6) 



The magnetic field is then derived from the Maxwell equation H = — x 
E/(//(a;)a;). The continuity conditions of the parallel components of the fields 
across the interface yield the following equations : 



Ex,i — — 0, (7) 
llE,,i + ^Ex,2^0. (8) 

We search a mode of the system which is a solution of the homogeneous 
problem. The system has a non-trivial solution if and only if 

At27i + /^i72 = 0. (9) 
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Taking into account equations (3,4), one obtains from (9) the surface wave 
dispersion relation for s-polarization: 



K 



2 _ u;^ AtiAt2[At2ei - Atie2] 
c2 - 



(10) 



For the particular case where ei 
simple form : 

^2 ^ /^i(a;)/^2(c^) 



/xi(a;) + //2(<^) 



62 = e, the dispersion relation takes the 



2.1.2 p-polarization (TM) 

For p-polarization, the electric field lies in the plane (y, z) and can be cast in 
the form: 



Ei(r,a;) 



E2(r,a;) 



^ ^ 

' ^ 
y^.,2y 



exp[i(K.R + 7i2;)], 



exp[i(K.R - 722;)]. 



(12) 



(13) 



The continuity of the tangential electric field yields 



(14) 



The Maxwell equation V.E = imposes a relation between the two compo- 
nents of the electric field 



(15) 



The continuity of the 2;-component of D yields: 

^\EzA — ^lEz-l- 



(16) 



Inserting (16) and (14) in (15) yields 
ei72 + £271 = 0. 



(17) 
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Taking into account equations (3,4), one obtains from (17) the surface wave 
dispersion relation for p-polarization: 

^ ^2— 2 ■ (18) 



For the particular case where iJ>i — IJ>2 — A*, the dispersion relation takes the 
simple form : 

= (19) 

ei(a;) + 62(0;) 



Remarks 

• When the media are non-magnetic, there are no surface waves in s-polarization. 
Indeed, the imaginary part of the z-components 7^ is always positive, so that 
7i + 72 cannot be zero. 

• At a material- vacuum interface (ei = /xi = 1), the dispersion relation reads 
in p-polarization 



K^\^^. (20) 
c\\e2{uj) + l 

It follows that the wave vector becomes very large for a frequency such that 
€{u;) + 1 = 0. 

The conditions (9) and (17) corresponds to the poles of the Fresnel reflection 
factors. To search these poles is an alternative and simple way to find the 
dispersion relation. This is particularly useful when searching the dispersion 
relation for multilayers system. 

For non-lossy media, one can find a real K corresponding to a real cu. This 
mode exists only if £2 < —1 in the case of an interface separating a vacuum 
from a material. 

In the presence of losses, the dispersion relation yields two equations but 
both frequency and wavector can be complex so that there are four pa- 
rameters. Two cases are of practical interest : i) a real frequency and a 
complex wavevector, ii) a complex frequency and a real wavector. These 
two choices leads to different shapes of the dispersion relation as discussed 
in [18,19,20,21]. The imaginary part of a; describes the finite lifetime of the 
mode due to losses. Conversely, for a given real cu, the imaginary part of K 
yields a finite propagation length along the interface. 

The dispersion relation (20) shows that for a real dielectric constant ei < — 1, 
K > uj/c. This mode cannot be excited by a plane wave whose wavevector is 
such that K < uj/c. In order to excite this mode, it is necessary to increase 
the wavevector. One can use a prism [14,22,23] or a grating [18]. A scatterer 
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can also generate a wave with the required wavevector. 
2.2 Dispersion relation 

In this subsection, we will consider two types of surface waves: surface-plasmon 
polaritons and surface-phonon polaritons. Surface-plasmon polaritons are ob- 
served at surfaces separating a dielectric from a medium with a gas of free 
electrons such as a metal or a doped semiconductor. The dielectric constant 
of the latter can be modelled by a Drude model : 

where ujp is the plasma frequency and F accounts for the losses. Using this 
model and neglecting the losses, we find that the resonance condition e{ui) + l ~ 
yields lo = uj^j \/2. For most metals, this frequency lies in the near UV so that 
these surface waves are more difficult to excite thermally. By contrast, surface- 
phonon polaritons can be excited thermally because they exist in the infrared. 
They have been studied through measurements of emission and reflectivity 
spectra by Zhizhin and Vinogradov [24]. Let us study the dispersion relation 
of surface-phonon polaritons at a vacuum/Silicon Carbide (SiC) interface. SiC 
is a non-magnetic material whose dielectric constant is well described by an 
oscillator model in the [2-22 /xm] range[25] : 

e{u^) = f 1 + ] , (22) 

with ujl = 969 cm'^ ujt = 793 cm-^ F = 4.76 cm"^ and = 6.7. The 
dispersion relation at a SiC/vacuum interface is represented in Fig. 2. This 
dispersion relation has been derived by assuming that the frequency ui is com- 
plex and the parallel wavevector K is real. This choice is well suited to analyse 
experimental measurements of spectra for fixed angles. The width of the res- 
onance peaks observed is related to the imaginary part of the frequency of 
the mode. We note that the curve is situated below the hght cone to — cK so 
that the surface wave is evanescent. We also observe a horizontal asymptote 
for uJasym = 1-784 10^"^ rad s^^ so that there is a peak in the density of electro- 
magnetic states. We will see in the next sections, that the existence of surface 
modes at a particular frequency plays a key role in many phenomena. 

In Fig. 3, we have shown the dispersion relation obtained when choosing a 
real frequency uo and a complex wavevector K. The real part of the complex 
wavevector is represented. It is seen that the shape of the dispersion relation 
is significantly changed. A backbending of the curve is observed. This type 
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Fig. 2. Dispersion relation for surface phonon-polariton at a SiC/Vacuum interface. 
The flat asymptote is situated at ujasym = 1-784 10^'' rad s~^. The slanting dashed 
line represents the light cone above which a wave is propagating and below which 
a wave is evanescent. 

of behaviour is observed experimentally when performing measurements at a 
fixed frequency and varying the angle. Observed resonances in reflection or 
emission experiments have an angular width which is related to the imaginary 
part of the complex wavevector. 

2.0x1 o" 



1.9 



~ 1.8 
'm 
T3 

to 



1.6 



1.5 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4x10° 

K(m" ) 

Fig. 3. Dispersion relation for surface phonon-polariton at a SiC/Vacuum interface. 
Real iv chosen to obtain a complex K. The real part of K is represented. The 
horizontal asymptote is situated at uJasym = 1-784 10^^ rad s~^. The slanting dashed 
line represents the light line above which a wave is propagating and below which a 
wave is evanescent - 

We have plotted in Fig-4 the surface wave decay length in the direction perpen- 
dicular to the interface versus the wavelength. From Eqs (1) and (2), it is seen 
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that the amphtude of the electromagnetic field decreases exponentially in the 
z direction with a decay length 8\ = 1/Jm(7i) in medium 1 and 62 = l//m(72) 
in medium 2. We note that the smallest penetration depth in SiC is obtained 




10.0 10.5 11.0 11.5 12.0x10 



Fig. 4. Surface wave decay length along the z-direction in medium 1 and 2 versus 
the wavelength for a SiC-vacuum interface. 

for the frequency ujasym- At this frequency, losses are very large. 

We study in Fig. 5 the surface wave propagation length along a SiC-vacuum in- 
terface. It is given by the inverse of the imaginary part of the parallel wavevec- 
tor L = l/Im{K). Around uJasym, L is minimum. It can be as large as several 
tens of wavelengths. 

It will be seen below that the existence of these surface modes is responsible 
for a long coherence time and a long coherence length of the electromagnetic 
field in the near field. They are essentially given by the lifetime and the decay 
length of the surface wave. 



1 -. 



0.1 - 




10.0 10.5 11.0 11.5 12.0x10' 

X (m) 

Fig. 5. Surface wave propagation length along the interface versus the wavelength 
A. 
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3 Fluctuation-dissipation theorem. Cross-spectral density. 

In this section, we introduce the tools and methods that are useful to derive 
the field radiated by a body in thermal equilibrium at temperature T both 
in the near field and in the far field. Whereas the phenomelogical theory of 
radiometry based on geometrical optics describes correctly the field emitted 
in the far field, it fails to predict the behaviour of the emitted radiation in 
the near field. Indeed, geometrical optics does not include evanescent waves. A 
new framework to describe thermal radiation is thus needed. Such a framework 
has been introduced by Rytov [26,27] and is known as fiuctuational electrody- 
namics. The key idea is that for any material in thermal equilibrium, charges 
such as electrons in metals, or ions in polar crystals undergo a random thermal 
motion. This generates fiuctuating currents which radiate an electromagnetic 
field. A body at temperature T is thus viewed as a medium with random 
currents that radiate the thermal field. The statistical properties of this field 
can be determined provided that i) the statistical properties of the random 
currents are known, ii) the radiation of a volume element below an interface is 
known. The first information is given by the Fluctuation-Dissipation Theorem 
(FDT), the second is given by the Green's tensor of the system. 

This approach is very similar to the Langcvin model for Brownian motion. 
Langevin [28] introduced a random force as a source for the dynamical equa- 
tions of the particles. This allows to derive the statistical properties of their 
random motion. An important feature of the model is that the random force 
is not arbitrary. Its correlation function is related to the losses of the system 
by the FDT. In the case of random electromagnetic fields, the dynamic equa- 
tions for the fields are Maxwell equations. We need to introduce an external 
random source to model the fiuctuations of the field. These external sources 
are random currents. The key issue now is to know the statistical properties 
of these random sources. They are given by the FDT. The remaining of this 
section introduces the technical tools needed to derive the radiated field. The 
first part is devoted to the statistical properties of the random currents given 
by the FDT, the second part deals with the FDT applied to the EM fields. 
The Green's tensor are given in appendix A. 



3.1 Cross-spectral density 

The spectral analysis of a signal is usually done using its Fourier transform. 

In the case of a stationary stochastic signal, the Fourier transform cannot 
be computed in the sense of a function because the integral is not square 
integrable. However, it is possible to compute the Fourier transform of the 
time-correlation function of the random signal. In what follows, we will be 
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interested in the space-time correlation function of the electromagnetic fields 
{Ek{r, t)Ei{r', t')). For a stationary field, this correlation function depends only 
on t — t' . Its Fourier transform £fc;(r,r', a;) is called cross-spectral density: 

oo 

£fcKr,r',a;)= j {Ek{y,t)Ei{y\t')) e'^^^'-^'Uit - t'). (23) 

— oo 

Note that for r = r', the above equation reduces to the Wiener- Khinchin 

theorem that relates the power spectral density of a random stationary signal 
to the Fourier transform of its time correlation function [29]. It is convenient 
to introduce a correlation function of the Fourier transforms using generalized 
functions : 

{Ek{v,u)El{v\J)) = 2T,5{uj-J)8u{r,r',uj). (24) 
3.2 Fluctuation- dissipation theorem for the current density 

The FDT derived by Callen and Welton [30] yields a general form of the 
symmetrized correlation function of a vector X(r,r', a;). Whereas for classi- 
cal quantities, the symmetrization does not change the results, it plays an 
important role in quantum mechanics for non-commuting observables. If the 
cross-spectral density of X is defined by A'fc/, we define the symmetrised cor- 
relation function of X by: 

'^fef^ = + ^ik]- (25) 

The symmetrised correlation function of the dipole moment of a particle in 
thermodynamic equilibrium with polarisability a defined by pj = e^aijEj can 
be written as [30,31]: 

{p,,{u;)p\{J))s = 2n5{uj-uj')v\!\uj) 

= hcoih{nLO /2kBT)Im[eoaki{u)]2T:5{uj - J), (26) 

where the brackets denote an ensemble average. In the preceding expression 
h is the reduced Planck constant and ks is the Boltzmann's constant. 

For a bulk in thermodynamic equilibrium at temperature T, the symmetrised 
correlation function of the polarization density can be written as [30,31]: 

(P,(r, uj)P:{v\ uj'))s = 2n6iuj - cu')vlp (r, r', to) 

^hcoth{hu;/2kBT)Im[eoeki{r, r', u)]2tt6{uj - uj'){27) 
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where the spatial dependence of the dielectric constant accounts for a possible 
non-locality. From this equation, we can easily derive the correlation function 
of the current density j = —iutP: 

{jk{r,u;)jl{r',u;'))s = Hu;'' coth{nu;/2kBT)Im[eoeki{r,r' ,u;)]2n5{u; - u;').{28) 



Note that if the medium is isotropic and local, the quantity Im[eki{r,r' ,00)] 
becomes Im[e{r,u;)]5ki5{r — r'). We also note that: 



^ coth{hu; /2kBT) = kv 



1 

+ 



2 e-K^{nw/kBT) - 1 



(29) 



is the mean energy of a harmonic oscillator in thermal equilibrium. In the 
following we will also use the compact notation 



for the mean energy of the harmonic oscillator without the zero point energy 

njjj/2. 



3.3 Fluctuation- dissipation for the fields. 



Another very useful application of the fluctuation-dissipation theorem yields 
a relation between the cross-spectral density of the fluctuating fields at equi- 
librium and the Green's tensor of the system [31]. The Green's tensor appears 
as the linear response coefficient relating the fields to their sources. Note that 
these quantities arc defined in classical electrodynamics. In what follows, we 
will use three different Green's tensors defined by: 



E 
H 



/<-> EE 
dh'G (r,r',a;)j(r',a;), (31) 

{v,u) = jd\'G (r,r',a;)j(r',a;), (32) 



and 



dh' G (r,r',u;)M(r',u;). (33) 
In this last equation G^i^ is the Green tensor relating the magnetic field to the 

^— > HE ■(r-¥ HH -(r-^ EE 

magnetization M. Note that both G and G are related to G through 
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the Maxwell equations [32]. 



G (r,r',a;) = -^Vr xG (r,r',a;), (34) 
G (r,r',a;) = -9^Vr X G (r, r', a;) x ^r'- (35) 

Here 

G^^(r,r',u;) x Vr' =^[Vr' G^^ {r,r' ,uj)] 

where the symbol T denotes a transposition of a tensor. Explicit forms of the 
Green's tensors are given in Appendix A. 

The cross-spectral correlation function of the electromagnetic field at equihb- 
rium or blackbody radiation then reads [31]: 

Sif\r,r',u) = i^ohw^coth[hw/2kBT]Im[G^i^{r,r',uj)]. (36) 

This last equation is the FDT for the electric field. The cross-spectral corre- 
lation function of the magnetic field obeys a similar relation 

nif (r, r', u) = eohuj^ coth[hu;/2kBT]Im[G^i^(r, r', u)] . (37) 

These relations yield the coherence properties of the equilibrium field provided 
that the Green's tensor of the system are known. Note that the Green's tensor 
is a classical object so that the vacuum fluctuations are already included in 
the above formalism. It is also important to note that the Green's tensor can 
be computed including the losses of the system. 

3.4 Relation between symmetrized correlation function and observables 

In what follows, we will use the FDT either with ©(a;, T) or with {huj/2) coih.[hw /2k bT] — 
Q{u!,T) + fuj/2. The former amounts to drop the vacuum energy hcu/2. This 
choice can be justifled using a heuristic argument [27]. When it comes to the 
derivation of fluxes, the vacuum energy cancels when taking the difference 
between emission and absorption [2 7]. Instead, when computing the Casimir 
force, one has first to compute the energy variation of the electromagnetic 
field in the space between two parallel plates. In that case, the vacuum fiuc- 
tuation energy hLo/2 plays a fundamental role and cannot be ignored so that 
it is kept in the calculation. This procedure may seem arbitrary. A more rig- 
orous approach to the choice of the relevant form of the FDT can be derived 
from quantum electrodynamics as discussed by Agarwal [31]. It can be shown 
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that when the process studied involves an absorption measurement, the rel- 
evant correlation function is the normally ordered correlation function [29]. 
If the measurement involves a quantum counter, then one needs to calculate 
the antinormally ordered correlation function [29]. The relevant forms of the 
FDT are given in the paper by Agarwal [31]. We show in Appendix B that 
one can end up with an effective cross-spectral density defined for positive 
frequencies only. This effective cross-spectral density depends on the type of 
measurement. The time-correlation function can be written as : 

{E,{v,t + T)E,{v\t)) = i?e[J^exp(^a;r)£,Y^(r,r',a;)], (38) 



where the effective cross-spectral density Elf^ {v, r', uj) that should be used for 
an absorption measurement is given by: 

4f)(r,r',^) = 4^/.o/m[Gf,^(r,r',a;)]e(a;,r), (39) 

and the cross-spectral density appropriate for a quantum-counter measure- 
ment is given by: 

£if\v, r', a;) = 4a;/.o/m[G'f,^(r, r', a;)][e(u;, T) + M- (40) 

Only the latter includes the energy of vacuum fluctuations. For the sake of 
comparison, we also report the symmetrised form apropriate for positive fre- 
quencies only : 

4?(r,r',a;) = 4a;/.o/m[Gf,(r, r', a;)][e(^, T) + ^]. (41) 

A simple rule can thus be used when starting with the usual symmetrised 
FDT as given by (36) in order to get the relevant correlation function for a 
process involving absorption: 1) restrict the spectrum to positive frequencies, 
2) multiply the spectrum by 2, 3) remove the energy fluctuation contribution, 
4) take the real part. 

3. 5 Fluctuational electrodynamics out of equilibrium 

In the previous section, we have given the form of the cross-spectral densities of 
the fields and current densities at equilibrium. However, it is possible to derive 
the fields radiated by a system out of equilibrium. The approach is based on 
the FDT for the current density. Assuming local thermal equilibrium, we can 
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derive the statistical properties of the currents. We can thus derive the fields 
radiated by a system with an inhomogeneous temperature field. Although 
the mean values of the fields are zero, their correlation are non-zero. Let us 
consider for instance the symmetrized cross-spectral correlation function of 
the electric field 

{E,ir,u)E*ir',u;'))s = 
(^l^luj'J dh,d'v^Gl^{v,v,,uj)Gl'^\v\v^,uj)jra^^^^ . (42) 

Using the FDT for the fluctuating currents (28), we obtain 

{E,{v,u)E;{v\J))s = ^ / dhjm[e{v,)] 
X [e[u;, r(rO] + ^/2]Gf^(r, n, uj)Gf^\v\ v,,J)2t:5{uj - J) . (43) 

With the help of the FDT, we have seen that it is possible to calculate all 
kind of cross-spectral spatial correlation functions involving the electric and 
the magnetic fields. With these functions, we are now able to calculate other 
quantities such as the energy density, the Poynting vector or the Maxwell stress 
tensor. In the case of thermal equilibrium situations, we will use the application 
of the FDT for the fields which give simpler expressions. Nevertheless, in non- 
equilibrium situation such as the study of heat transfer between materials held 
at different temperature, these expressions are no longer valid. It is however 
still possible to use the fluctuation-dissipation theorem for the currents by 
assuming local thermal equilibrium. It will thus be possible to derive the fluxes 
for non-equilibrium situations. 



4 Electromagnetic energy density and Local Density Of States 
(LDOS) 

In this section, we will study how the electromagnetic energy density is mod- 
ifled by the presence of material media. We shall first examine the amount 
of electromagnetic energy emitted by a half-space at temperature T. It will 
be shown that the density of energy is dramatically different in the near field 
and in the far field when surface waves are excited. The second point that we 
address is the general problem of the definition of the local density of electro- 
magnetic states. Whereas it is possible to derive the density of electromagnetic 
states for a non-lossy system by searching the eigenmodes, the lossy case is 
more difficult. An alternative approach was introduced by Agarwal [33] based 
on the FDT. Using this approach, we will discuss the role of surface waves. 
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4-1 Density of emitted electromagnetic energy 

Let us return to the configuration described in Fig.l. We calculate the elec- 
tromagnetic energy density in the vacuum above a material (medium 2) at 
temperature T due to the emission of this material. We do not take into ac- 
count the energy incident on the medium. In order to retrieve the density of 
energy at equilibrium, we should include the energy incident on the surface. 
The density of energy in vacuum reads (See Jackson [34] p. 242) 

(t/) = |(|E(r,t)r} + f (|H(r,t)p) = /^KU^,^), (44) 



where we have introduced a spectral density of energy utot- The details of the 
derivation are given in [7,35]. The basic procedure amounts to derive the field 
radiated by the random currents in the lower half-space. Adding the electric 
and magnetic contributions, the total electromagnetic energy above a medium 
at temperature T in a vacuum at T = K is 




where the Fresnel reflection factors are given in the appendix. 
4-2 Discussion 

In order to illustrate this discussion, we study the density of electromagnetic 
energy above some specific materials. Let us first consider a material support- 
ing surface waves in the infrared such as SiC. In Fig. 6, we plot the energy 
density Utot{z,uj) versus the frequency at different distances of a semi-infinite 
interface of SiC. The semi-infinite medium is at temperature T = 300 K 
whereas the vacuum is at T = K. Note that at T = 300 K, Wien's law gives 
a peak wavelength for thermal radiation Xwien — 10//m. In the far field, i.e. 
for distances d larger than AvKjen/27r, the energy density spectrum resembles 
that of a blackbody. The difference with a Planck spectrum comes from the 
fact that SiC is a very refiecting material around A = lO/xm or a; = 1.7 10^^ 
rad s~^. Thus, its emissivity is small in this frequency interval. This property 
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100 200 300 40C500x10^^ 
(0 (Hz) 

Fig. 6. Total electromagnetic energy density above a plane interface separating SiC 
at T = 300 K from vacuum at T = K. 

is easily recovered from the electromagnetic energy due to propagating waves 
only (First term in (45)) 

/• //O 1 _ Ir* P -I- 1 _ Ir^ |2 

= «°(^,T) / ^LJlj^L^LJlJiL^ (46) 



where we have used 2'KKdK = cos{9)dfl, 9 is the angle between the emission 
direction and the nomal of the surface. The integral is performed over a half- 
space and 



T) = ^ ^ = ^e(c., T) (47) 

TT^c^ exp ^ — 1 c-^ 
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is the electromagnetic energy density in a cavity at thermal equilibrium T. 
In the far field, the evanescent waves do not contribute to the energy density 
because of the exponential decay (e~^^™('>'i)^). We note that if medium 2 is 
totally absorbing (r^f* = 0); the energy density due to propagating waves is 
half the energy calculated in a vacuum at thermal equilibrium. This is not 
surprising since we are computing only the emitted part of the radiation. In 
the case of equilibrium radiation, there is also the contribution of the radiation 
coming from the upper half space. We note that in (46), the emissivity appears 
to be (1 - |r^2(^'^)P + 1 - ki2(^^,^^)n/2. It is thus the half sum of the 
energy transmission factors for both polarizations. Thermal emission by a 
half-space can be viewed as a transmission process of a blackbody radiation 
in the material medium through an interface. This point of view yields insight 
in Kirchhoff's law. Indeed, the equality between emissivity and absorptivity 
appears to be a consequence of the equahty of the transmission factor when 
interchanging source and detector. 




100 200 300 400 500x10^^ 

M (Hz) 

Fig. 7. Electromagnetic energy density above a plane interface separating glass at 
T = 300 K from vacuum at T = K. 

At a distance z = 3/xm which is slightly larger than Avyien/2vr, the energy 
density spectrum changes drastically and a strong peak emerges. At 100 nm 
from the interface, one observes that the thermal emission is almost monochro- 
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matic around lo — 1.78 lO^^rad s~^. At this frequency the energy density has 
increased by more than four orders of magnitude. The peak corresponds to 
the excitation of a surface wave. This distance is in agreement with the decay 
length of the surface waves as discussed in Section 2. At distances much smaller 
than the wavelength, we enter a regime that we call extreme near field. The 
leading contribution comes from the very large wavevectors K in the energy 
density integral. At large it can be shown that 71 ~ iK ^ so that lm{r\,^ 
tends to zero and /m(rf2) tends to its electrostatic limit [e^ — l)/(e2 + !)• We 
finally obtain a very simple asymptotic form of the energy density 



showing that the energy density will diverge at the frequencies where e\ ~ 
— 1. These are the frequencies where the dispersion relation has horizontal 
asymptotes. For lossy materials, a resonance will occur at the frequency where 
Re{e2) = —1 provided that absorption (i.e Jm[e2]) is not too large. Let us 
mention that the resonance of the reflection factor in the electrostatic limit 
has been observed experimentally by Hillenbrand et a/. [3]. A spectrum of the 
light scattered by a tip very close to a surface shows a peak for the resonance 
frequency. This peak is due to the field emitted by the image of the tip. Indeed, 
its amplitude is proportional to the refiection factor (e2 — 1)/ (£2 + 1) and is 
therefore resonant when Re{e2) = —1. 

In order to prove that this surprising behaviour is not specific to SiC, we 
plot in Fig. 7, the energy density spectrum above a flat interface of glass. This 
material is amorphous so that the optical phonons are poorly defined. Here 
again, the energy density spectrum resembles that of a blackbody in the far 
field whereas peaks emerge while approaching the surface. The strongest peak 
is at cij = 9.2 10^^ rad s^^ and the weakest one at c<j = 2.2 10"'^'^ rad s~^. These 
frequencies are solution of Re[e2{uj)] = —1. Yet, the corresponding surface 
waves have a very short propagation length. 

Not all materials supporting surface waves exhibit strong peaks in their near- 
field thermal energy density spectrum. Indeed, as it can be seen in (48), a 
peak is exhibited if the frequency where €2(00) approaches -1 corresponds to a 
frequency range where Q{uj,T) is not too small. For example, metals exhibit 
surface-plasmon polariton in the UV or visible range where Q{uj,T) is expo- 
nentially small at ambient temperature. Thus metals do not exhibit strong 
peak in their thermal energy density spectrum in the near field. 
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4-3 Local Density Of States 



The density of states (DOS) is a fundamental quantity from which many 
macroscopic quantities can be derived. In statistical physics, the DOS allows 
to calculate the partition function of a system from which all the macroscopic 
properties follow. The local density of states (LDOS) is useful to study a non- 
uniform system. The local density of electronic states is widely used in solid 
state physics. It has been shown [36] for instance that a scanning tunneling 
microscope images the electronic LDOS. The local character of the LDOS 
clearly describes the spatial distribution of electrons in the solid. A similar 
spatial dependence is also relevant for electromagnetic waves. Whereas the 
intensity is uniform in a vacuum in equilibrium, this is not the case in a 
waveguide or above an interface. The distribution of the energy is no longer 
uniform. Whereas LDOS is well defined for electrons in solid state physics [37], 
its electromagnetic counterpart is not very well defined in the literature. As 
compared to electronic systems, two differences must be taken into account : 
the vectorial nature of the fields and the existence of losses. 

In electrodynamics, the LDOS is used in different contexts. From Fermi's 
golden rule, it is known that the DOS determines the radiation rate. It can be 
shown that the lifetime of an atom with an electric dipole along a unit vector 
u is inversely proportional to /m[u ■ G^^ ■ u]. This is often rcfcrcd to as the 
LDOS. To avoid confusion, we shall refer to this quantity as the projected 
LDOS. It is the relevant quantity that one needs to study when designing a 
microcavity or a photonic crystal to tailor emission properties. Yet, note that 
only those states that can be coupled to the dipole are taken into account. 
Thus it is not a good definition if one is interested in the total energy of the 
system. Such a quantity is required when computing dispersion forces [38,39] 
or shear forces [40] for instance. Those forces depend on the energy stored 
in all available modes. In a vacuum, the LDOS can be shown to be given by 

the imaginary part of the trace of the Green's tensor G ■ This seems to 
be a straightforward extension of the scalar result which is proportional to 
the imaginary part of the Green function. The vacuum form is thus usually 
assumed to be valid for any other situations [41,42,43]. In what follows, we 
will summarize a recent analysis of the LDOS [44] that follows the original 
approach by Agarwal [33]. It will be seen that the LDOS is not given by the 

imaginary part of the trace of the Green's tensor G . It will also appear that 
surface waves dominate the LDOS close to an interface. 

We consider a system in thermal equilibrium at temperature T. In a vacuum, 
one can define the electromagnetic energy U{u) by the product of the DOS 
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p{u!) by the mean energy of each state at temperature T 



(49) 



Wc can now introduce [7,44] a local density of states by using as a starting 
point the local density of electromagnetic energy U (r, cu) at a given point r in 
space, and at a given angular frequency cu. This can be written by definition 
of the LDOS p(r, u;) as 

U{r,u;) = p(r,a;) ^ - 

exp[nu!/kBT) — 1 



(50) 



The density of electromagnetic energy is the sum of the density of electric 
energy and of the density of magnetic energy given in (44). In equilibrium, 
it can be calculated using the system Green's function and the fluctuation- 
dissipation theorem. We start from the electric and magnetic field correlation 
functions for a stationary system 



{Ek{r,t)Ei{r',t'))^Re 

{Hk{r,t)Hi{r',t')) = Re 



^4r(r,r',a;)e-(-') 



/|^<'(r.r'..)e- 

.0 



■iu){t-t') 



(51) 

(52) 



with t = t' . The cross-spectral density for normally ordered fields is given by 
(39). It follows that the energy per unit volume can be cast in the form [33,44]: 



[exp{hu}/kBT) - 1] TTC^ 



ImTr 



^EE ^HH 

G (r,r,a;)-FG (r,r,a;) 



.(53) 



A comparison of Eqs. (50) and (53) shows that the LDOS is the sum of an 
electric contribution and a magnetic contribution p^: 



p{r,uj) 



-ImTr 



TTC^ 



G (r,r,u;) + G (r,r,a;) 



= p^(r,a;) + p^(r,a;)(54) 



In what follows, we shall discuss a few examples to illustrate the modification 
of the LDOS. It will be seen that in some cases, the LDOS is accurately given 
by the trace of the electric Green's dyadic but it can also be very different. 
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4-4 Electromagnetic LDOS in simple geometries 



4-4-^ Vacuum 

In the vacuum, the Green's tensors G and G obey the same equation 
and have the same boundary conditions. Therefore, their contribution to the 
electromagnetic energy density are equah 

4^ EE i—^HH (jj ^ 

Im[G (r, r, oj)] = Im[G (r, r, a;)] = — I . (55) 

DTTC 



The LDOS is thus obtained by multiplying the electric field contribution by 
2. After taking the trace, the usual result for a vacuum is retrieved 

p^(r, uj) = py{uj) = (56) 



As expected, we note that the LDOS is homogeneous and isotropic. 



4.4-2 Plane interface 

We now consider a plane interface separating a vacuum (medium 1 in the upper 
half-space) from a semi-infinite material (medium 2, in the lower half-space) 
characterised by its complex dielectric constant 62 (c^). The material is assumed 
to be homogeneous, linear, isotropic and non-magnetic. The expression of the 
LDOS at a given frequency and at a given height z above the interface in 
vacuum is obtained by inserting the expressions of the electric and magnetic 
Green's tensors for this geometry [45] into equation (54). Note that in the 
presence of an interface, the magnetic and electric Green's tensors are no longer 
the same. Indeed, the boundary counditions at the interface are different for 
the electric and magnetic fields. 

Let us consider some specific examples for real materials like metals and di- 
electrics. We first calculate p{z^ uo) for aluminum at different heights. Alu- 
minum is a metal whose dielectric constant is well described by a Drude model 
for near-UV, visible and near-IR frequencies 



u){uj -\- iV) 



with Up = 1.747 10^6 rad.s"^ and F = 7.596 10^^ rad.s"^ We plotted in Fig.8 
the LDOS p(r, a;) in the near UV-near IR frequency domain at four different 
heights. We first note that the LDOS increases drastically when the distance to 
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Fig. 8. LDOS versus frequency at different heights above a semi-infinite sample of 
aluminum. From [44] 

the material is reduced. As discussed in the previous paragraph, at larger dis- 
tances from the material, one retrieves the vacuum density of states. Note that 
at a given distance, it is always possible to find a sufficiently high frequency for 
which the corresponding wavelength is small compared to the distance so that 
a far-field situation is retrieved. This is clearly seen when looking at the curve 
for z = 1/im which coincides with the vacuum LDOS. When the distance to 
the material is decreased, additional modes are present: these are the evanes- 
cent modes which are confined close to the interface and which cannot be seen 
in the far field. Moreover, aluminum exhibits a resonance around to = cUp/ a/2. 
Below this frequency, the material supports surface-plasmon polaritons so that 
these additional modes are seen in the near field. This produces an increase of 
the LDOS close to the interface. The enhancement is particularly important at 
the resonant frequency which corresponds to Re[e2{uj)] = —1. This behaviour 
is analogous to that previously described in the electromagnetic energy density 
paragraph for a glass surface supporting surface-phonon polaritons. Also note 
that in the low frequency regime, the LDOS increases. Finally, Fig. 8 shows 
that it is possible to have a LDOS smaller than that of vacuum at some par- 
ticular distances and frequencies. Fig. 9 shows the propagating and evanescent 
waves contributions to the LDOS above an aluminum sample at a distance of 
10 nm. The propagating contribution is very similar to that of the vacuum 
LDOS. As expected, the evanescent contribution dominates at low frequency 
and around the surface-plasmon polariton resonance, where pure near-field 
contributions dominate. 

Let us now turn to the comparison of p{z,Lj) with the usual definition often 
encountered in the literature[41,42,43], which corresponds to p^{z, u). We plot 
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Fig. 9. Density of states contributions due to the propagating and evanescent waves 
compared to the total density of states and the vacuum density of states. These 
quantities are calculated above an aluminum sample at a distance of 10 nm. Prom 
[44] 




Fig. 10. LDOS at a distance z = 10 nm above a semi-infinite aluminum sample. 
Comparison with p^{z,uj) and p^{z,uj). From [44]. 



in Fig. 10, p, and p^ above an aluminum surface at a distance z = 10 nm. 
In this figure, it is possible to identify three different domains for the LDOS 
behaviour. We note again that in the far-field situation (corresponding here to 
high frequencies i.e. X/27i <C z), the LDOS reduces to the vacuum situation. 
In this case p{z, u) = 2p^{z, u) = 2p^{z, u). Around the resonance, the LDOS 
is dominated by the electric contribution p^ . Conversely, at low frequencies, 
p^{z,uj) dominates. Thus, Fig. 10 shows that we have to be very careful when 
using the approximation p{z,uj) = p^{z,uj). Above aluminum and at a dis- 
tance z = 10 nm, it is only valid in a narrow range between u = 10^^ rad s~^ 
and uj = 1.5 x 10^^ rad s~^, i.e. around the frequency where the surface wave 
exists. 
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4.4-3 Asymptotic form of the LDOS in the near-field 

In order to get more physical insight, we have calculated the asymptotic LDOS 
behaviour in the three regimes mentioned above. As we have already seen, the 
far- field regime (A/27r -C z) corresponds to the vacuum case. To study the 
near-field situation (A/27r ^ z), wc focus on the evanescent contribution due 
to the large wavevectors K as suggested by the results in Fig.9. In this (quasi- 
static) limit, the Fresnel refiection factors reduce to 



Asymptotically, the expressions of p^{z,u!) and p^{z,u!) are[44]: 



Pv 



e2 + l|2 4A;g^3' 



p"{z,uj)^py 



+ 



8koz 2|e2 + iP/co^ 



(60) 
(61) 



At a distance z = 10 nm above an aluminum surface, these asymptotic expres- 
sions matches almost perfectly with the evanescent contributions {K > ko) of 
p^ and p^ . These expressions also show that for a given frequency, one can 
always find a distance z to the interface below which the dominant contribu- 
tion to the LDOS will be the one due to the imaginary part of the electric- 
field Green function that varies like {koz)~^. But for aluminum at a distance 
z = 10 nm, this is not the case for all frequencies. As we mentioned before, 
this is only true around the resonance. For example, at low frequencies, and 
for z — 10 nm, the LDOS is actually dominated by pye2 / (IGkoz) . 



4 -4 -4 Spatial oscillations of the LDOS 

Let us now focus on the LDOS variations at a given frequency u versus the 
distance z to the interface. There arc essentially three regimes. First, for dis- 
tances much larger than the wavelength, the LDOS is given by the vacuum 
expression p^. The second regime is observed close to the interface where os- 
cillations are observed. Indeed, at a given frequency, each incident plane wave 
on the interface can interfere with its reflected counterpart. This generates an 
interference pattern with a fringe spacing that depends on the angle and the 
frequency. Upon adding the contributions of all the plane waves over angles, 
the oscillating structure disappears except close to the interface. This leads to 
oscillations around distances on the order of the wavelength. This phenomenon 
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is the electromagnetic analog of Friedel oscillations which can be observed in 
the electronic density of states near interfaces [12,47]. For a highly reflecting 
material, the real part of the reflection coefficients are negative so that the 
LDOS decreases while approaching the surface. These two regimes are clearly 
observed for aluminum in Fig. 11. The third regime is observed at small dis- 
tances as seen in Fig. 11. It is due to the contribution of surface waves. Its 
behaviour is thus dependent on the frequency. Let us first consider the partic- 
ular case of the frequency corresponding to the asymptote of the dispersion 
relation. It is seen that the evanescent contribution dominates and ultimately 
the LDOS always increases as l/z^, following the behaviour found in (60). 
This is the usual quasi-static contribution that is always found at short dis- 
tance [35] . At a frequency slightly lower than the resonance frequency, surface 
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Fig. 11. LDOS versus the distance z from an aluminum- vacuum interface at the 
aluminum resonant frequency. Prom [44] 

waves are still excited on the surface. These modes increase the LDOS accord- 
ing to an exponential law, a behaviour which was already found for thermally 
emitted fields [8,35]. 

Conclusion about the LDOS 

The main results of this section can be summarized as follows. The LDOS 
of the electromagnetic field can be unambiguously and properly defined from 
the local density of electromagnetic energy in a vacuum above a sample at 
temperature T in equilibrium. The LDOS can still be written as a function of 
the electric-field Green's tensor only, but it is in general not proportional to 
the trace of its imaginary part. An additional term proportional to the trace 
of the imaginary part of the magnetic-field Green's tensor is present in the far 
field and at low frequencies. At short distance from the surface of a material 
supporting surface modes (plasmon or phonon-polaritons), the LDOS has a 
resonance at frequencies such that Re[e{uj)] = —1. Close to this resonance, 
the approximation p{z,uj) = p^{z,u) holds. If it is possible to measure the 
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near-field thermal emission spectrum of a material, the local density of states 
can be retrieved [44]. 



5 Coherence properties of planar thermal sources in the near-field 

In this section, we examine the second order coherence properties of the fields 
due to thermal excitation in the presence of surface waves. We have shown that 
the density of energy is completely dominated by the contribution of surface 
waves in the near field. We shall see that they are also responsible for a deep 
modification of coherence properties. In what follows, we restrict ourselves to 
second-order coherence properties. 

5.1 Spatial coherence in the near field 

The spatial coherence of the electromagnetic field is characterized by its cross- 
spectral density Ekiij^i '^)- Roughly speaking, we study the correlation func- 
tion of the electromagnetic field at two different points for a particular fre- 
quency. For a system in thermal equilibrium, this quantity is readily given by 
the fiuctuation-dissipation theorem (36). The spatial dependence is thus in- 
cluded in the spatial dependence of the imaginary part of the Green's tensor. 
For the particular case of a vacuum, one finds the properties of the blackbody 
radiation. The Green's tensor is given in appendix A. When taking its imagi- 
nary part, it is seen that the spatial dependence of the cross-spectral density is 
proportional to sin(A;r)/A;r. The typical length scale of the correlation function 
for the blackbody radiation is thus the wavelength A = 27i/k. 

It is more interesting to analyse what happens in the presence of an interface. 
It turns out that the coherence length may be either much larger or much 
smaller than the wavelength. This problem has been studied in [8,35] for a 
slightly different case. The authors considered the coherence properties of the 
field emitted by a sofid. The difference with the above result is that in equi- 
librium, one has to consider two contributions : i) the blackbody radiation 
illuminating the surface and reflected by the surface, ii) the radiation emitted 
by the surface. In what follows, we focus only on the emitted part of the ra- 
diation so that we do not consider the equilibrium situation. The correlation 
is given by (43). This equation is valid inasmuch as the temperature T can 
be defined everywhere in a half-space. It requires a local thermal equilibrium. 
We have represented in Fig. 12 the cross-spectral density of the electric field 
for different metallic surfaces at a given distance z = 0.05A to the interface. 
It is seen that the correlation oscillates and has an exponentially decaying en- 
veloppe. The decay length is much larger han the wavelength indicating that 
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Fig. 12. Cross spectral density function £'a;a;(ri, r2, a;) (denoted Wxx in the label of 
the figure) of the thermally emitted x component of the electric field versus p/X 
where p = |ri — for different metals (a) and for SiC at different wavelengths (b). 
The long range correlation is due to surface plasmon polaritons for metals and to 
surface phonon polaritons for SiC. From [8] 

the fields are coherent over large distances. This surprising phenomenon is due 
to the excitation of surface waves along the interface. The physical mechanism 
is based on the fact that a small volume element contains random currents 
that excite a surface wave. This surface wave propagates along the interface 
over distances larger than the wavelength. It follows that different points may 
be illuminated by the same random source so that they are correlated. Accord- 
ingly, one does not expect any correlation between the s-polarizcd field since 
no surface wave exists for s-polarization. If one uses a material with a real 
part of the dielectric constant larger than —1, no surface wave can propagate 
so that no correlation should be observed. We have shown in Fig. 12 the case 
of tungsten in the visible that does not support surface waves. It is seen that 
the coherence length is smaller than a wavelength so that the radiation field 
appears to be more incoherent than blackbody radiation. It has been shown 
in [35] that in the short distance regime, the coherence length is given by z. 

A similar behaviour is observed for SiC, a polar material that supports surface- 
phonon polaritons in a frequency band. Within this band, at a wavelength 
11.36//m the correlation is seen to be a long-range correlation whereas the 
correlation decays very rapidly for a wavelength (9.1//m) that is not in the 
band where surface waves exist. 

Let us discuss in simple terms the physical origin of these unusual coherence 
properties. The long-range coherence is unexpected because the fluctuating 
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currents are delta correlated as shown by the FDT. This is the reason why 
the fields are usually assumed to be delta-correlated in space. However, the 
fluctuating currents excite induced currents in the material. In the case of 
a metal, a surface plasmon can be excited. In the case of a polar crystal, 
a surface-phonon polariton can be excited. Both surface waves are extended 
modes along the surface. The induced currents associated with these extended 
modes are therefore coherent over large distances. More precisely, the coher- 
ence length is expected to be given by the decay length of these surface modes. 
This has been confirmed by a detailed asymptotic analysis in [35] . 

The other surprising property shown in Fig. 12 is that the coherence length 
defined as the FWMH of the cross-spectral density can be smaller than the 
wavelength. In other words, a source can be more spatially incoherent than 
the blackbody radiation. The key idea is that close to an interface, the field 
contains evanescent waves so that features smaller than the wavelength can 
exist. This is not the case in a vacuum so that the field has a minimum 
coherence length. Since the amplitude of evanescent waves of large wave vector 
K decays as exp{—Kz), it is clear that the distance z appears as a cutoff 
frequency. This explains that the coherence length increases as z in the near- 
field regime. 



5.2 Temporal coherence in the near field 



The temporal coherence of the electromagnetic field is characterized by the 
same point time-correlation function of the electromagnetic field : 

{Ek{r,t + T)Ei{r,t)). (62) 



This correlation function is a measurement of the memory of the random 
field. It is useful to introduce a typical decay time tcoh of the correlation func- 
tion called coherence time. A Michelson interferometer with aligned mirrors 
performs a measurement of the correlation function. Indeed, the interference 
term of the signal can be written as E{r,t + T)E{r,t) where r is the flight 
time corresponding to the optical path length difference Sgpt between the two 
paths 5 opt = CT. If the path- length difference is larger than the longitudinal 
coherence length ctcoh, no interferences can be observed. 

The temporal coherence of the EM field is related to its power spectral density. 
This is clearly seen by using the Wiener-Khinchin theorem [29,48] which shows 
that the power spectral density is the Fourier transform of the correlation 
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function. Alternatively, we can start from (36). It follows that 



{Ek{r,t + T)Ei{r,t))^Re 



I A^ioooeiu;,T)Irn[G^f^iv,v,u;)]e'^^^ 

.0 



.(63) 



Let us first consider the temporal coherence of the field in a vacuum. The 
imaginary part of the Green's tensor does not diverge and yields zero for 
non diagonal terms and cu/Gttc for diagonal terms. It follows that the time 
correlation function of the blackbody radiation is given by : 



{Ek{r,t + T)Ei{r,t))^Ski Re 



1 A^,uQ{u,T)^e^-^— 



(64) 



Since the integrand has a large spectral width, it appears that the coherence 
time is on the order of the peak radiation period. 

If we now consider the case of an interface, we know that the spectrum can 
be very different in the near field. We have seen previously that the contribu- 
tion of the surface wave modifies dramatically the density of electromagnetic 
energy. In particular, we have seen that the density of energy becomes quasi- 
monochromatic which suggests a large coherence time. More specifically, in 
the extreme near field, we have seen in (59) that the Green's function has a 
resonant denominator e + 1. Close to the resonance where i?e[e(a;o)] = —1, we 
can expand the dielectric constant as : 



e a; 



-1 + «e"(tUo) + i}^ - t^o) 



de 
do;' 



(65) 



where we have used the notation e — e' + ie" . The denominator e + 1 can be 
cast in the form: 



e{uj) + l^[ — \[uj-ujQ + iVl (66) 




where V = e"(c<Jo)/-^- It is seen that the Green's dyadic has a pole at the 
frequency corresponding to the asymptote of the dispersion relation of the 
surface wave. Its contribution to the integral (64) yields an exponential de- 
cay of the form exp[iu!ot — Ft]. It follows that in the extreme near field, the 
thermally emitted field is partially temporally coherent with a coherence time 
given by F"^. The origin of the temporal coherence of the electromagnetic field 
can thus be assigned to the very large density of states due to the asymptote 
of the surface wave. It follows that whereas the plane interface of a hot metal- 
lic surface is a temporally incoherent source for an atom located in the far 
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field, it is a partially temporally coherent source for an atom located within a 
nanometric distance from the interface. 



5.3 Polarization coherence in the near field 

We have seen that the excitation of surface waves by delta-correlated currents 
produces both a spatial and a temporal coherence. The correlation of the field 
is characterized by the decay time and the decay length of the surface wave 
that propagates along the interface. There is a further interesting coherence 
effect that has been studied by Setala et al. [49]. The basic idea is that surface 
waves are p-polarized waves. As shown in section 2, a surface wave that prop- 
agates along the y axis has an electric field with two components along the y 
and z axis. It follows that the y and z components of the field are correlated. 

The study of correlation of the fields in the near field of the source cannot be 
performed using the standard formalism. Indeed, when dealing with a beam, 
it is usually possible to work with the two components of the field perpen- 
dicular to the propagation axis. This is no longer possible in the near field 
of a thermal source. A generalization of the degree of polarization has been 
introduced in [49]. It was found that the degree of polarization varies as a 
function of the distance to the interface. Like the coherence time, it increases 
when approaching the interface from the far field because the surface wave 
creates a correlation. When reaching the very near-field regime, the degree of 
polarization decays and tends to 1/4 for all the materials. 



6 Spatially peirtially coherent thermal sources in the far field 

6. 1 Introduction 

In this section, we will discuss the possibility of designing a source that is 
spatially partially coherent. In simple terms, a spatially partially coherent 
source is a source that radiates a field which has a narrow angular aperture at 
a given wavelength. The typical examples of coherent sources are lasers and 
antennas. These sources have well defined emission angular lobes. In what 
follows, we will show that a narrow angular emission lobe is a signature of the 
spatial coherence of the field in the plane of the source. 

We first introduce an analogue of the Wiener Khinchin theorem (WKT) (63) 
to analyse the spatial coherence of the field. In simple terms, WKT states 
that a quasi monochromatic source with bandwidth Ai/ has a coherence time 
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roughly equal to 1/Az/. Similarly, a quasi parallel beam with a spatial fre- 
quency bandwidth Ak^ has a transverse coherence length 27r/AA';^.. A formal 
proof is based on two properties : the relationship between the cross-spectral 
density of the field in the plane of the source z — and the power spectral 
density of the field, the relationship between the power radiated in far field 
and the power spectral density. For a translationally invariant system, the 
Fourier transform of the field docs not converge in the sense of a function. 
Yet, one can define a field equal to the random field in a square of area A and 
null outside. We can now define the Fourier transform of the field in the plane 
2; = as : 

E^(r||, a;) = j Ea(K, u) exp(iK • r,,) (67) 



It can be shown [48] that the WKT yields a relation between the cross-spectral 
density and the power spectral density of the spectrum of the field : 

{Ei{^ , uj)E:{^ + rJi , ^)) = I ^ jim 1 (e,,^(K, ^)£;,%(K, u)) exp(zK • rj) (68) 



This relation implies that if the spectrum has a bandwidth smaller than 27r/A, 
the coherence length is larger than A. The second step is to show that the band- 
width of the spectrum power density <^E;^^(K, a;)E;*^(K, a;)^ /A is given by the 
emission pattern in the far field. Indeed, the field can be written everywhere 
as [50,29] : 

EA(r, ^) = / E^(K, exp(iK • r,, + i-fz), (69) 



where k = (K, 7) and 7 is given by ^'^ — kQ—K.'^. This field can be evaluated 
asymptotically in the far field using the stationary phase approximation [29]. 
It can be cast in the form : 

E^(r, u) = ^^^-^ E^ [K = — r II , a; j , (70) 



where is a constant. The power dP flowing through an element of surface 
dS = r^dfl is given by the flux of the Poynting vector. In far held, the Poynting 
vector has locally a plane wave structure so that its time averaged amphtude 
is given by eoc|Ep/2. 



dP='-^\K\' 



E. 



A 



dfl 



(71) 



where f is the unit vector r/|r|. This relation completes the discussion of the 
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link between the directivity of the emitted field and the coherence of the field 
in the plane of the source. It is clear that a directional source implies a narrow 
spectrum and therefore a large correlation length. 



6.2 Design of coherent thermal sources 

We have seen in the previous section that a source which supports a surface 
wave is partially spatially coherent along the surface. However, because these 
waves cannot propagate in a vacuum, the coherence remains confined in the 
vicinity of the surface. The question that we address now is whether it is 
possible to export the near-field coherence to the far field. In essence, that 
amounts to couple the surface waves to the propagating waves. This can be 
done in several ways. A practical way is to rule a grating on the surface. 
The grating can then diffract the surface wave. By properly choosing the 
period of the grating, it is possible to control the angle of propagation of the 
diffracted light. This was first observed in [51,52,53] for a very deep grating 
ruled on a doped silicon surface. Such a material supports surface-plasmon 
polaritons in the infrared. A more effective source was realized using a gold 
grating by Kreiter et al. [54]. Heinzel et al. [55] have also reahzed a source 
using surface plasmons on tungsten in the near infrared. A source based on 
the use of surface-phonon polaritons on SiC was reported by Lc Gall et al. 
[18]. The first discussion of the spatial coherence of these sources was reported 
in [10]. An extended discussion of these properties has been recently reported 
[56] where transverse coherence lengths have been derived from experimental 
measurements of angular widths of emission peaks. Angular peaks as narrow 
as 1° can be obtained. 

In order to have an efficient source of light, it is necessary to make sure that 
the coupling of the surface wave into a propagating wave is optimum. This 
can be accomplished by designing a surface with total absorption. From Kirch- 
hoff's law, it follows that if absorption is unity, then emissivity is also unity. 
Another condition must be satisfied : the emission wavelength should lie in 
a region where Planck's function takes large values. The optimum choice of 
the wavelength thus depends on the temperature of the source. Fig. 13 shows 
the angular emission pattern of a SiC grating. It is clearly seen that the an- 
gular aperture is very narrow indicating a large coherence length [56]. Let us 
emphasize that the coherence is not due to to the grating but only to the 
surface wave. The role of the grating is merely to couple the surface wave into 
a propagating wave. 

Different schemes have been proposed to produce partially coherent thermal 
sources. They are based on a filtering of the emission pattern in order to re- 
duce the angular width of the emission pattern. The scheme that we have 
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Fig. 13. Experimental angular emission of a SiC grating at two different wavelengths. 
The angular pattern has the characteristic shape of an antenna. It demonstrates the 
spatial coherence of the thermal source. Measurements are taken at 800 K. Prom 
[56] 

described so far is based on the use of a grating that couples a surface wave 
to propagating wave can be viewed as a device that increases the absorptiv- 
ity or emissivity to 1 for a narrow set of angles. It can thus be viewed as a 
filtering process. A different type of filter can be designed using multilayers 
systems. Several authors [57,58] have reported a narrow angular pattern emis- 
sion obtained by interferences between several layers. This mechanism leads 
to angular widths on the order of 10°. 

6. 3 Engineering radiative properties of surfaces 

For many applications, it is desirable to modify radiative properties of sur- 
faces. An introduction to radiative properties of surfaces can be found in a 
review paper by Z. Zhang [59]. Roughness has often been used to increase the 
emissivity. An analysis of the different mechanisms involved can be found in 
[60] . Further references on scattering by rough surfaces can be found in several 
reviews and monographs [50,61,62,63,64,65]. Microstructures can be used to 
design efficient selective absorbers and sources. The decay of reficctivity of a 
shallow rough surface due to the excitation of surface plasmons is addressed 
in [66]. Hava and coworkers have examined silicon microstructured surfaces 
[67,68,69,70]. Sai et al. [71] have designed silicon microstructured surfaces for 
thermophotovoltaics applications. Marquier et a/. [72] have studied the effect 
of surface plasmon on highly doped silicon showing that the peak absorption 
frequency can be tuned by varying the doping. Kusunoki et al. [73] have re- 
ported emissivity measurements on tantalum surfaces with two-dimensional 
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periodic structures. They observed peaks of emission due to the excitation of 
surface plasmons. Pralle et al. [74] have designed selective infrared emitters 
using periodic structures on sihcon wafers coated with gold. 




Fig. 14. Angular emissivity measurements of a SiC grating at the surface phonon 
polariton frequency. The emissivity is above 0.9 and quasi-isotropic in the plane of 
incidence. The sample temperature is 800K. Prom [56]. 

An interesting application of surface waves to the enhancement of light emis- 
sion has been demonstrated by Marquier et al. [56,72]. The idea is to emit 
light in all directions at a given frequency. To this aim, it is necessary to be 
able to couple light in all directions. When using a grating of period d, the 
emission angle 6 is given by ^ sin(^^) = K{uj) + p'^ where p is an integer. 
The key idea is to work at the resonance frequency taking advantage that the 
dispersion relation is fiat. At that particular frequency, surface plasmons exist 
for a very large set of wavevectors K. It follows that one can design a source 
emitting in all directions. This has been demonstrated for a SiC surface on 
which a grating was ruled. Fig. 14 represents the angular emission pattern. It is 
seen that the emissivity is close to its maximum value 1 and almost isotropic. 
The mechanism of emission can be again viewed as a two steps process. First, 
each volume element is equivalent to a random dipole that can either emit a 
photon or excite a surface wave. The latter is a very efficient process so that 
usually, most of the desexcitation goes into surface waves and eventually into 
heat. The second step is the introduction of a grating that converts the surface 
wave into a propagating wave by diffraction. Thus, the excitation of a surface 
wave which usually tends to decrease the emission of light becomes a factor 
that enhances the emission of light. A similar mechanism has been proposed 
recently to use the surface waves in order to enhance the emission of light 
by quantum wells placed close to a metallic layer [75,76]. The idea is that the 
high density of states due to the surface plasmon enhances the emission. If the 
conversion of a surface wave into a propagating wave is efficient, the process 
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enhances the emission. 



7 Radiative heat transfer in the near field 

We have shown previously that the density of electromagnetic energy increases 
in the near field due to the contribution of surface waves. We now address the 
question of heat transfer between bodies separated by distances smaller than 
the wavelength. In that case, it is expected to observe contributions of the 
surface waves to the radative heat transfer. This topic has already a long his- 
tory. Anomalous radiative heat transfer was observed in the sixties. Cravalho 
et al. [77], Bochm and Tien [78] studied that problem and took into account 
waves that are propagative in the materials and evanescent in the vacuum 
gap. Yet, this correction does not take into account all the evanescent waves. 
Waves that are evanescent on both sides of the interface (i.e. surface modes) 
were not taken into account in these early works. The first correct derivation 
of the fiux between two plane parallel plates was reported by Polder and Van 
Hove [79] in 1971. Their method allowed them to analyse the contributions of 
different polarizations and to compute separately the contributions of evanes- 
cent and propagating waves. Similar works were reported later on by different 
authors [80,81,82,83,84,85]. It is found in all these works that the flux diverges 
as the distance decreases. This point was considered to be unphysical in [86,87] 
where the increase in the density of states due to surface waves was not in- 
cluded. It is pointed out in [88,89] that when keeping the contribution of all 
the wavevectors, the flux diverges if the temperature difference is assumed to 
be kept to a constant value. This amounts to say that the thermal resistance 
goes to zero which docs not violate any physical law. Only a few experiments 
[90,91,92] were reported on measurements of heat transfer due to near-field 
radiation. 

Heat transfer between a plane and a small particle was first discussed by Doro- 
feyev [93] and later by Pendry [82] and Volokitin and Pcrsson [83]. Practical re- 
sults were derived for a metal using a Drude model and making the additional 
assumption that |e| >> 1. Mulct et al. [94] pointed out the resonant contribu- 
tion of surface waves to the heat transfer. It was shown that the heat transfer 
is quasimonochromatic at the frequency of the optical phonons for a polar 
crystal given by Re[e{uj)\ + 1 = 0. A similar effect is observed between metal- 
lic parallel surfaces [85] and for doped semiconductors [72]. For metals, this 
resonance does not play a significant role because the plasma frequency is in 
the UV domain so that the Bose- Einstein factor takes low values at usual tem- 
peratures. It was later suggested [10,95] that a quasi-monochromatic enhanced 
heat transfer could be used to increase the efficiency of thermophotovoltaics 
devices by matching the energy of the emitted photons with the absorption 
band gap of the photovoltaics cell. That might reduce the loss of excess energy 
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of ultraviolet photons. The heat transfer between two small particles has been 
studied in [83]. It has been shown that the dipole-dipole interaction yields a 
large contribution to the heat transfer whereas the contribution of the photon 
emission and absorption process is neghgible. This near-field heat transfer be- 
tween nanoparticles is analogous to the energy transfer between molecules due 
to the dipole-dipole coupling known as Forster transfer [96]. It may also be 
resonant for surface plasmon resonances. In what follows, we shall derive ex- 
plicitly the heat transfer between nanoparticles and the emission by a surface 
plane. 



7.1 Radiative power exchanged between two spherical nanoparticles 



We now calculate the heat transfer between two spherical nanoparticles held 
at different temperature and separated in the vacuum. Such a calculation was 
first reported by Volokitin and Persson [83] . Let us consider two nanoparticles 
whose dielectric constant are ei and t2 and whose temperatures are Ti and T2. 
We first calculate the power dissipated in particle 2 by the electromagnetic 
field induced by particle 1 using the dipolar appoximation: 

Pi_2H = eo|/m(Q;2)|Ei„e(r2,a;)|', (72) 



where Y2 denotes the position of the particule 2 and a2 is the polarisability of 
a sphere of radius a [97]: 

a2 = 47ra^^^^. (73) 
62 + 2 



The field incident on the particle 2 created by the thermal fiuctuating dipole 
of particle 1 located at ri at temperature Ti is given by: 

Ej„c(r2,cu) = /xoCt;^G(r2,ri,tc;) • p (74) 



where G(r2! ri, cj) is the vacuum Green's tensor given in appendix A. To pro- 
ceed, we need the correlation function of the dipole given by the fiuctuation- 
dissipation theorem whose symmetrised form is given by (27). We finally ob- 
tain the heat exchange between two spherical nanoparticles at temperature Ti 
and T2. 

^1-2 = ^ 16 i®(^' Ti) - e>(^, T2)\ . (75) 
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Let us note the spatial depcndancc of the heat transfer. This dependence 
is typical of the dipolc-dipolc interaction. It is actually a van der Waals type 
interactions that can be interpreted in the following way: fluctuations (thermal 
or not) distort the charge distribution of a nanoparticle producing a fluctuating 
dipole. This fluctuating dipole induces in turn an electromagnetic fleld on the 
other nanoparticle initiating a second dipole. This dipole interaction causes 
both an energy transfer and a momentum transfer or force. For molecules 
this energy transfer is known as Forster transfer and the force is called van 
der Waals force. We flnd that nanoparticles follow a similar behaviour with a 
resonance at the surface polariton resonance. Indeed, in the case of spherical 
particle with radius a, the polarisability has a resonance when the dielectric 
constant approaches —2 provided that the imaginary part of the dielectric 
constant is not too high. The particle resonances appear in the visible part of 
the spectrum for metals and in the infrared for polar materials. 



7.2 Thermal emitted flux by a planar interface 



In this section, we analyse the emission by a plane interface. Let us consider 
the situation of a planar interface {z = 0) separating a dielectric {z < 0) at 
temperature T from a vacuum > 0). We shall derive the flux emitted from 
the interface. 



7.2.1 Classical theory of radiation. 

In the classical theory of radiation, the power d'^Q emitted by an elementary 
opaque surface dS at temperature T, in a solid angle dQ around a direction u 
making an angle 6 with the normal to the surface (Fig. 15), whose monochro- 
matic emissivity is e'(^) is 

d''Q{u;,e)^eMlUT)dndS (76) 



where 



A-K^c^ Q^I^bT _i 

is the blackbody speciflc intensity. The power dQ emitted by the elementary 
surface is thus 

dQ{u) = y e'^{Q)ll{T) cosedndS ^q-ndS (77) 
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Fig. 15. Solid angle dQ around a direction u making an angle 9 with the normal of 
an elementary surface dS. 

where the integral is performed over a half-space. We have introduced the 
radiative vector defined by 

qH= Js'MiT)udn (78) 

In the case of a blackbody, the integration over the angles is straightforward. 
The monochromatic heat flux is thus equal to 7rL° (T) and the total heat flux 
is equal to (tT^ where a is the Stefan constant 

a = = 5.67 lO-'Wm-^K-^ 

7.2.2 Fluctuational electrodynamics method 

We will follow Polder and van Hove [79] to show that the phenomenologi- 
cal form of the emitted flux can be derived in the framework of fluctuational 
electrodynamics. In the following sections, we use this approach to derive the 
expressions valid in the near field. Let us consider the situation of a planar 
interface {z — 0) separating a dielectric {z < 0) at temperature T from a 
vacuum > 0). The flux emitted by the half-space is given by the Poynting 
vector S = E X H. In the case of monochromatic quantities, the time aver- 
age Poynting vector reads S(r,cj) = |i?e (E(r, c<j) x H*(r,u;)). This quantity 
can be derived using the fluctuation-dissipation theorem. Thus, the electric 
and magnetic Green's tensor are needed. In this geometry, these tensors are 
given by (A. 3) and (A.4). The Poynting vector reduces to its z component 
{S^{r,uj)) = {l/2)Re[(E^Hl - EyH*)]. In order to obtain the Poynting vec- 
tor, one calculates quantities like ^£'i(r, a;)if*(r, o;)^. Using (31,32) 

{e,{v,u;)H*{v,u;)) = (^/ioc / ^.^/(r, ^^^^^^(r, r")i,(r')ir(r")rfVdV) (79) 
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Using the effective FDT for the currents (28) defined for positive frequencies 
only, the preceding equation reduces to 



E,{r,uj)H*{r,uj) 



TTC^ 



/6"(r')Gf/(r,r')Gf/*(r,r')dV (80) 



Using the Green functions expressions (A. 3) and (A. 4) and the identities e" = 
in the upper half space and that ^I^LtP' !<? = 2i?e(72)/m(72), the Poynting 
vector can be cast in the form 



{S^{r,uj))= „ — Re 



X 



\j.s |2 I |2 

r2il "T r2il 
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(81) 



Let us note that only 71 may not be real in the preceding expression. In fact, 
as the upper half space is vacuum (medium 1), 71 can only be real (K < cu/c) 
or pure imaginary (K > uj/c). Therefore, the contribution of evanescent waves 
vanishes in the radiative flux expression. Using the identities (82-85), 



it:e(6*7)-i^e(7) '^' 7, 



Im{e*^) — Im{'j) 



12I71I 



= i?e(7i)(l - IrliH - 2/m(7i)/m(r|0 



I72I 



(82) 
(83) 
(84) 
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Re{e;^2)Kf= i?e(ehi)(l-|r^ir)-2/m(ehi)MrfO .(85) 



p |2^ 



one finally obtains: 



LU/C 



KdK 1 - |r 



12 



(86) 



As already mentioned, only propagating waves {K < u/c) contribute to this 
expression. This is not surprising because no waves come from the positive 
z direction. Moreover, there is a revolution symmetry around the z axis. In- 
troducing dVt the elementary solid angle, we have the relation KdK/kQ — 
dfl cos 9 /{2Tr). The Poynting vector is then given by 



bT-1 / 



cosOdQ 1- |r^2p + l 



\rP 12 
\'12\ 



Q. = 2-K 



27r 



(87) 
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In the case of a blackbody, i.e. a body for which the reflection factors are 
null, the usual expression of the radiative flux 7r/°(T) is recovered. When the 
dielectric situated below the interface does not behave as a blackbody, the flux 
takes the usual form 

qiu) = {S,{r,u;)) ^ j du e'Ml^{T) cosOdQ (88) 



where we have identified the emissivity e'^{9) = (1 — |r*2p + 1 — Ir^aH/S. 
In presence of a single interface, we note that the radiation emitted is not 
different from the usual one. The near field does not play any role in this 
situation. 



7.3 Heat transfer between two semi-infinite polar materials. Interference ef- 
fects 



We now focus on the heat transfer between two semi-infinite half-spaces sep- 
arated by a vacuum gap and whose temperature Ti and T2 are uniform. The 
main changes that occur at small distance is the fact that evanescent waves 
can contribute to the heat transfer through tunneling. 

We summarize in the next sections the results. Detailed derivations can be 
found in [85] for instance. The radiative transfer is characterised by the radia- 
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Fig. 16. Two semi-infinite half-spaces separated by a vacuum (distance d) 
tive flux. In the phenomenological theory, this flux is given by 

cosOd^l I duj- — 

P'lU!P'2u! 



(89) 
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where the e^^ are the directional monochromatic emissivities and the p'^ are 
the directional monochromatic reflectivities. 



Using a fluctuational electrodynamics model, the flux can be written as the 
sum of two terms qiiS) — (f^'^iiS) + (f'"°''^{oj). 

The first term g^'"°^(a;) is the propagating waves contribution. 
= ^ I dwd^ cose 

q=s,p ^ 



(1 














1 - 







CiTi) - I'^{T2)] (90) 



Let us note that 1 — l^^sf P and 1 — \r^f\'^ are the transmission energy coeffi- 
cients between media 1 and 3 and 2 and 3 for the s or p-polarization. These 
coefficients can be identified as an emissivity in the same way that it has 
been defined for a single interface. Let us remark that this expression for the 
propagating waves contribution to the radiative flux between two semi-inflnite 
media is very close to the usual one. Only the denominators are different be- 
cause interferences are not taken into account in the phenomenological model. 
Nevertheless, if one considers a frequency interval small in comparison with the 
frequency but larger than c/ d, the variation of e*'^*'^ with uj is much faster than 
the Presnel reflection coefficient variations. The integration over this interval 
yields an average value of |1 — r\^r\2G^^^^'^\'^ which is exactly 1 — |'"3inr32p. 
Matching the reflectivity with the Fresnel reflection energy coefficient, one can 
then identify this expression for the radiative flux with the classical one. 

7.3.1 Tunneling of evanescent waves 

The second term g^^""(c<;) is the contribution of the evanescent waves. It reads: 



q=s,p •', ^0 

u>/c 

Contrary to the single interface case, this contribution does not vanish because 
of the existence of both upward and downward evanescent waves in the space 
between the two media[98]. When the distance reduces, this term is more 
and more important due to the presence of the exponential e''^^M'y3)d_ \Yhen 
the material involved are supporting surface waves, the imaginary part of the 
reflection coefficient in p-polarisation becomes important around the resonant 
frequency, when the dielectric constant approaches -1. If the two media are 
sufficiently close to allow the interaction between the exponentially decaying 
surface waves bound to each interface, a transfer occurs due to the tunneling 
of evanescent waves. 



CiTi)-Cm] (91) 
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Fig. 17. Radiative heat transfer coefficient versus distance for semi-infinite media of 
temperature T = 300 K. From [85] 
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Fig. 18. Monochromatic heat transfer coefficient for a distance d = 10 nm and a 
temperature T = 300 K. From [85] 

Let us define a radiative heat transfer coefficient as the hmit of the ratio of 
the radiative flux on the temperature difference between the two media when 
this temperature tends to zero: 



hm 



q{uj) 



(Ti-T2)^0 Ti - T2 



(92) 



In Fig. 17, h^{(jj) is represented versus the distance between two semi-infinite 
media of glass or SiC. For a distance larger than the thermal radiation wave- 
length given by the Wien law, i.e. for c? > 10 /xm, the transfer does not depend 
on the distance. We are then in the classical regime where the transfer occurs 
through the radiation of propagating waves. At shorter distances, the transfer 
increases as For a distance of 10 nm, the radiative heat transfer coeffi- 

cient has increased by 4 orders of magnitude compared to its far field value. 
If we now focus on the spectral dependance of the heat transfer coefficient 
at a 10 nm distance (Fig. 18), we note that the heat transfer is important for 
the frequencies corresponding to resonant surface waves. The heat transfer 
is therefore practically monochromatic in the near field. We can also expand 
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asymptotically the radiative heat transfer coefficient for short distances: 

,p, , 1 /m(ei)/m(e2) , ( V e^^/^^^^^ 



This expression yields the l/d? depcndance of the transfer coefficient and 
its strong frequency dependance. Indeed, when the dielectric constant ap- 
proaches -1, the radiative heat transfer coefficient exhibits a peak as well as 
the Fresnel reflection factor. This is the signature of the presence of a surface 
wave. The validity of the dependence has been questionned in [86] on 
the basis that an infinite flux is not physical. As a matter of fact, the flux is 
infinite if one assumes that the temperature difference is kept constant. This 
problem is analogous to the problem of an electron flux or intensity that goes 
to infinity if the resistance goes to zero at fixed voltage. This raises another 
question: in the case of heat transfer, the resistance across a vacuum gap on 
the order of 10 nm usually remains much larger than the bulk conduction 
resistance of solids over distances on the order of 100 nm. It is thus safe to 
assume that the temperature is uniform in the bulk over a skin depth so that 
the calculation is valid. 



7-4 Calculation of the heat transfer between a dielectric sphere and a half- 
space 



7.4-1 Introduction 

We calculate in this part the radiative power exchanged between a small spher- 
ical particle and a semi-infinite medium. To this aim, we first calculate the 
power absorbed by the dielectric sphere placed above a heated half-space. We 
then calculate the power dissipated by the half-space situated below a heated 
sphere from reciprocity [94]. The geometry of the problem is presented in 
Fig. 19 : the upper medium 2; > is vacuum. A particle (P) of radius a and 
dielectric constant ep{uj) = ep(cj) + iep{uj) is held at temperature Tp. The 
lower medium is filled by a homogeneous, isotropic material (bulk) of dielec- 
tric constant eB{w) = ^b{'^) +^%('^) ^-^d held at temperature Tb- The center 
of the particle is at a distance d above the interface. 



7.4-2 Power absorbed by the bulk- Near-field limit 

A calculation following the procedures already introduced (see e.g. [83,94]), 
yields the power absorbed by the particle when illuminated by the field radi- 
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Fig. 19. Geometry of the system 

ated by a half-space: 




If we consider the fluctuating currents inside the particle that radiates into 
the bulk and dissipates, one can calculate by the same formalism the power 
locally dissipated per unit volume at a point r inside the bulk. It reads 

P^^r^'lr,^) = -^/m[eBH]7m[aH]e(u;,rp))^|G„^(r,rp,u;)|2 (95) 

^ ^ n,m 

From the expression of the one-interface Green's tensor, it is possible to expand 
asymptotically the expression of the power absorbed by the particle. This 
quantity behaves as and reads 

From this expression, we see that there is an enhancement of the power ab- 
sorbed if the denominators vanish or approach zero. We have seen that it is 
the case if the material support resonant surface waves so that the dielectric 
constant of the material can take negative values. We study in the next section 
the case of SiC. 
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7.4-3 Example of SiC 



As it has been said in the first part of the article, SiC is a polar material that 
can be described by an oscillator model (22). In Fig.20, we plot P^J^^(a;) for 
a spherical particle held at temperature Tp = 300 K, of radius a = 5 nm at 
different distances above the surface. We note that the figure displays two re- 
markable peaks at frequency uji ^ 1.756 x 10^"^ rad.s^"*^ and 002 ~ 1.787 x 10^^ 
rad.s^^ These two peaks correspond to the resonances of the system. The 
first one corresponds to a frequency where ep{uj) approaches -2 : a volume 
phonon polariton is excited in the particle inducing a large electric dipole 
and a large dissipation. The second one is related to the resonant surface 
wave corresponding to a large increase of the electromagnetic LDOS. Thus, 
here too, the radiative heat transfer in the near field can be considered as 
monochromatic. The electromagnetic waves associated with the resonant sur- 




Fig. 20. Mean power radiated by the bulk (at Tg = 300 K) and absorbed by the 
particle (of radius a=5 nm) versus frequency, (a) d=20 nm; (b) d=50 nm; (c) d=100 
nm. The insert (log-log scale) shows the spectrum of the absorbed power between 
10^^ rad s~^ and 10^^ rad s"^; (e) d=20 nm; (f) c?=lmm. Prom [94] 

face waves are evanescent. The energy transfer, which finds its origin in the 
presence of these waves, is important because the particle lies in the region 
(up to many micrometers) where the evanescent field is large, so that there 
is an efficient coupling between the field and the particle. In the far-field, 
evanescent waves are negligible and usual results are retrieved. Fig. 21 shows 
the integrated power absorbed by the same particle versus the distance d. The 
near-field radiative heat transfer increases as 1/d^ (as it was suggested by the 
asymptotic behavior) and is larger at small distances by several orders of mag- 
nitude than the far field one. This enhancement comes from the contribution 
of evanescent waves. Reciprocity requires that the same enhanced radiative 
heat transfer appears when the particle illuminates the surface. This situation 
may help to understand the radiative heat exchange between a nano-tip (like 
those used in near-field microscopy) and a sample. To answer this question, we 
calculated from (95) the total power (integrated over frequencies) dissipated 
per unit volume for different points in the sample. Fig. 22 displays a map, in 
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d in m 



Fig. 21. Total power radiated by the bulk (at Tb = 300 K) and absorbed by the 
particle (of radius a = 5 nm) versus distance. Prom [94] 

log-scale, of the dissipation rate in the case of a lOnm-diameter sphere of SiC at 
Tp — 300 K situated 100 nm above a sample of SiC. It is seen that the energy 
is dissipated on a scale comparable to the tip-sample distance. The dissipation 
per unit volume decreases very fast (1/r^) with the distance r between the 
source and the point of the sample where the dissipation is considered (the 
isocontour labeled with a '6' corresponds to the points where the dissipation 
per unit volume is 10^ W m~^. The amount of energy locally deposited is as 
large as 100 MW m~^. We note that the dependance of the heat deposited 
follows the same regime as in the two particles exchange in the vacuum. This 
is not surprising. Here also the phenomenon is due to an interaction between 
induced dipoles. This phenomenon is also at the origin of the force between 
macroscopic bodies at nanometric distances which is the subject of the next 
section. 




Fig. 22. Deposited power per unit volume inside the bulk. The particle has a radius 
a = 5 nm and is held at temperature Tp = 300 K. Prom [94] 
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8 Role of surface electromagnetic waves on the Casimir force 

8. 1 Introduction 

After having considered the energy exchange due to the interaction in the 
near field between surface waves, it is natural to wonder what happens in 
terms of momentum exchange in the near field when two semi-infinite bodies 
are approached face to face. This situation is actually well known since 1948, 
when Casimir [99,100] first showed the existence of an attracting force between 
two parallel perfect conductors. A large body of litterature has been devoted 
to this effect and several reviews are available [101,102,103,104,105,106]. The 
seminal paper of Lifshitz [107] occupies a special place because it was the 
first calculation of this force by means of the FDT for the currents. Agarwal 
[108] reported a similar calculation using the FDT for the fields. There has 
been an increasing interest in the Casimir force since it has been shown that 
this force could be measured with high accuracy [109,110,111,112,113,114] and 
that it should be considered in the design of micro electromechanical systems 
(MEMS) [115,116]. Various correction to this force have been studied such as 
finite conductivity [117] or temperature corrections [118,119]. In what follows, 
we discuss the role of surface waves in the Casimir force. We will show that 
they play a key role in the short distance regime. 

To get a simple picture of Casimir force, we recall that a system of two plane 
parallel refiecting planes is a waveguide. The number of electromagnetic modes 
in the waveguide is discrete and depends on the thickness of the waveguide. 
From quantum electrodynamics, it is known that each mode with frequency uj 
has a minimum energy fvjj/2 refered to as vacuum fiuctuations. If the thickness 
decreases, the number of modes decreases so that the electromagnetic energy 
decreases. Hence, the existence of vacuum fluctuations entails that there is an 
attractive force between the two plates. This phenomenon is clearly a macro- 
scopic manifestation of the electromagnetic energy of the vacuum which is a 
pure quantum effect. Yet, its computation amounts to count the number of 
electromagnetic modes available and this is a pure classical problem. As has 
been discussed in the previous sections, the local density of electromagnetic 
states is completely dominated by the existence of surface modes. It follows 
that the role of surface waves is essential in the physics of the Casimir force 
in the short distance regime. 

In classical electrodynamics, the momentum transfer is given by the Maxwell 
stress tensor Tij [34]. The flelds can be derived using the FDT [120]. Using 
the system Green tensors for the electric and the magnetic fleld, it is possible 
to obtain this quantity. In the case of two semi-infinite bodies separated along 
the z axis by a vacuum gap, the momentum fiux reduces to the zz component 
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Tzz of the Maxwell stress tensor given by 



T 

J- 7. 



fo 
2 



\Err, 



l-Ej/l — \Ez\ 



+ 



l-f^xl^ + \iiy^ ~ \Hz 



(97) 



One has to substract the infinite contribution to the force in the absence 
of bodies [121]. The force can be attractive or repulsive depending on the 
materials properties [122,123]. One obtains an attractive force in the case of 
dielectrics [124] and a force that might be repulsive in some configurations 
implying magnetic materials [125,126]. 



8.2 Spectrum of the force 



Lifshitz [107] obtained a force per unit area given by 



oo oo 

du rdu 



^2.7 2.^^^'^) ^''^ 

where v = {v? — 1)^/^ (Imv < 0), and is the Fresnel refiection coefficient 
for a plane wave with polarization and wavevector K = parallel to the 
vacuum-medium interface. We use he convention that an attractive force cor- 
responds to F < 0. The force appears as the contribution of elementary plane 
waves which angular frequency is u and which wavevector parallel to the inter- 
face is ULu/c. In his paper, Lifshitz used a deformation contour in the complex 
plane of frequencies to obtain a final formula where the summation over the 
wavevector is replaced by an integral over v and the summation over the fre- 
quencies is replaced by an integral over the imaginary frequencies u = i^. 
This approach has the advantage to replace the oscillating exponentials with 
smooth real functions that make the integral easy to integrate numerically. 
Nevertheless, by doing such a deformation contour, one is losing the spectral 
information contained in the expression (99). What we are going to show in 
the following is that the main contribution to the force in the near field is 
coming from the coupled polaritons of both interfaces. Therefore, we will see 
that there is a complete analogy in the interpretation of the momentum ex- 
change with the interpretation of the energy exchange in terms of interaction 
of surface-polaritons. 

Let us study then the force spectrum in the case of two real materials. In 
the case of SiC (Fig. 23), the force is dominated by the UV and IR contribu- 
tions. Actually, due to the presence of the u;^ in the expression (99), the UV 
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contribution is much more important than the IR one. 
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Fig. 23. Contributions of s and p-polarized, propagating and evanescent modes to 
the force spectrum (as given by (99) after integration over the wavevector u). 
Distance d = 10 nm. Material: SiC, dielectric function with two resonances. The 
angular frequencies of the corresponding surface resonances are 1.78 x 10^^ s"^ in 
the IR and 2.45 x 10^^ g-i ^he UV [46]. Prom [120] 



7 2 

X 



-4 







prop, s 

- - ■ prop. p 

evan. s 

evan. p 







'iV lo" lo" lo" lo" lo" lo' 



angular frequency (rad.s ) 



Fig. 24. Contributions of s and p-polarized, propagating and evanescent modes to 
the force spectrum ( (99) integrated on the wavevector u). Distance d = 10 nm. 
Material: Aluminum, described by tabulated optical data [46]. From [120] 

Fig. 24 represents the force spectrum between two aluminum half-spaces. It 

is seen that all the contributions are coming from frequencies very close to 
plasmon-polariton resonances. From expression (99), it is easy to separate the 
polarisation contributions, the propagative contribution (0 < m < 1) and the 
evanescent one (m > 1). We see in Figs. 23, 24 that the peaks exist only in p- 
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polarisation and for evanescent waves. This fact is an additional argument in 
favor of an interpretation of the force being caused by the interaction between 
polaritons which are not propagating and only appear in p-polarisation. An 
interpretation of the phenomenon is the following. At frequencies close to the 
system resonances i.e. close to the surface waves frequencies, the density of 
electromagnetic states increases. Thus, the amount of momentum carried in 
the gap increases too ^ . This interpretation is only valid at distances of the 
order of the wavelength of the surface wave. At larger distances, the force is 
dominated by the propagative contribution as the original perfect conductors 
situation first pointed out by Casimir. We note also that the contribution to 
the force at such small distances has a different sign whether the frequency is 
lower or higher than the resonant frequency. However, the total Casimir force 
between dielectrics is always positive and smaller than the one for perfect 
conductors. 

8.3 Binding and antibinding resonances 

The role of SPP can be further analysed by studying the variation of the 
integrand of the force F{u, uj) in the plane {u, uj). Close to the resonant surface- 
wave frequency Wswt the contributions to the force come from evanescent 
waves. We therefore limit our study to the case m > 1 and close to a; ~ 
cusw- In Fig. 25a, we plot the integrand F{u,lj) for two aluminum half-spaces 
separated by a distance d = 10 nm. Two branches mainly contribute, the 
higher frequency branch yielding a repulsive contribution whereas the lower 
one gives an attractive contribution. These two branches actually follow the 
two-interfaces system dispersion relation given by 

1 - rle-^"^-" = (100) 

The infiuence of the dispersion relation on the force is illustrated in Fig. 25b. In 
this figure, the quantity 1/|1 — r^e~'^~^'^\'^ is plotted in the {u,uj) plane. Upon 
comparison between Fig. 25b and Fig. 25a, it is clearly seen that the main 
contribution to the force can be attributed to the SPP. In addition, we observe 
on Fig. 25b a dark line which corresponds to minima of 1/|1 — TpC^^c^'^p. The 
minima can be attributed to very large values of the refiection factor of a plane 
interface r^. Thus, the dark line is the dispersion relation of a single SPP on 
a flat interface. In Fig. 26, the integrand is plotted for two aluminum half- 
spaces separated by a distance d = 100 nm : the two branches tend to merge 
with the fiat interface dispersion relation. One can thus propose the following 
interpretation : when the surfaces approach each other, the overlapping of 

^ As for energy transfer case, momentum is carried by evanescent waves if and only 
if an upward and a downward wave are present. 
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Fig. 25. (a) Wave-vector resolved spectrum of the Casimir Force (99) in the {u, lu) 
plane between two aluminum half-spaces separated by a distance of 10 nm. The 
frequency of the flat asymptote corresponds to the peaks of the force spectrum 
(Fig. 24). Light (dark) areas : attractive (repulsive) force, (b) Resonant denominator 
|1/1 — rpe~^~'"^p in the {u,lu) plane, the grayscale giving the logarithm to base 10. 
The dispersion relation of the coupled surface resonnance correspond to light areas; 
dark area : dispersion for a single interface (20) . The dielectric function is extracted 
from tabulated data [46] . The inset sketches the magnetic field of the coupled surface 
resonances (antisymmetric and symmetric combinations). From [120] 

the two SPP leads to two coupled modes and to a splitting of the polaritons 
frequencies [127,128]. The frequency splitting can be found from the solutions 
of (100) which are implicitely defined by 

rp(M,u;) =±e'"^'^/^ 

The signs correspond to either symmetric or antisymmetric mode functions 
(for the magnetic field), as sketched in the inset of Fig. 26. The symmet- 
ric (antisymmetric) branch corresponds to a lower (higher) resonance fre- 
quency, respectively, similar to molecular orbitals and tunneling doublets. 
These branches contribute with opposite signs to the Casimir force, due to 
the following identity: 

2r2(cu,M)e-2'^^'^/'^ 
1 -r2(a;,it) q-'^^-^^-Ic ~ 
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Fig. 26. Same as Fig.25 but for d = 100 nm. From [120] 
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(101) 



where the first (second) term is peaked at the symmetric (antisymmetric) 
cavity mode. The symmetry of the resonance mode function hence determines 
the attractive or repulsive character of its contribution to the Casimir force. 



8.4 Analytical formulation of the short- distance limit 



Using a simple Lorentz-Drude model for the dielectric function, 

= , , ^^e?!^. (102) 



we can derive a simple analytic form of the force in the near-field limit that 
accounts for the modification of the local density of states due to the surface 
waves and that takes into account absorption effects. This model describes 
either dielectric or metals depending on the value of cuq. The corresponding 



plasma frequency is — ciJq). With this convention, the large u asymptote 

of the SPP dispersion (20) occurs at ~ ^2. As the distance d reduces to a 
quantity small compared to the wavelength, the evanescent contribution to 
the force comes from higher and higher parrallel wave vector. When ^ 1, 
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the reflection coefficient in p-polarisation can be approximated by: 



> 1 : 



rp{uj,u) 



- 



— i'juj — uj^ 



(103) 



It is seen tliat tlie reflection factor fias a complex pole. Prom (100), we thus 
get the following dispersion relation for the (anti)symmetric surface-plasmon 
resonances, neglecting for the moment the damping coefficient 7: 



(104) 



where we have used w ~ m for m 3> 1. For large u, we solve by iteration and 
find that u!±^ft. As announced above, the symmetric mode thus occurs at a 
lower resonance frequency. 

To derive an analytical estimate for the Casimir force, we retain in (99) 
only the contribution of p-polarized, evanescent waves. The final result[120] 
contains two contributions: 



F = 
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47rfi 



(105) 



where 
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(106) 
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(107) 



When (2; — > 1), we get the following asymptotics 
a{z) ^ 0.1388 - 0.32(1 - z) + 0.4(1 - zf 
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X (1-^)2 (109) 



with ({3) ~ 1.202. For a metal, the parameter z takes the value 1 so that we 
have a(l) ^ 0.1388 and Lig ^ 1.202. We compare (105) in Fig.27 to the full 
integral (99) for the case of aluminum : the asymptotic estimation turns out 
to be quite accurate for distances d < O.lXspp where Xspp is the wavelength 
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Distance (m) 

Fig. 27. Comparison of different expressions for the Casimir force between alu- 
minum surfaces. We plot the ratio F{d) / Fcas{d), where Fcas{d) = Hctt'^ / {2A0d'^) 
is the Casimir Force for perfect mirrors. Solid line: numerical integration of (99), 
using tabulated data [46]. Short-dashed line with circles : same, with a model di- 
electric function of Drude form (102) with ujq = 0, Q = 1.66 x 10^^ s~^, and 
^/Vt = 0.036. These parameters have been derived by fitting the value of the reflec- 
tivity. Long-dashed line: short distance asymptotics (105) with the same values for 
o^o, ^ and 7. From [120] 

of the SPP with the largest frequency. In the case of aluminum, the first order 
correction in 7/!] is 2.5% of the zeroth order value of the force. The plot also 
shows that for numerical integration, the tabulated data[46] and the Lorentz- 
Drude model (102) with parameters fitted around the surface resonance give 
very close results over a large range of distances. This is another indication 
that the short range Casimir force is dominated by a narrow frequency range. 



8. 5 Friction forces 

The Casimir force that has been treated in the previous section is not the only 
one that occurs when two bodies are approached one to each other. A particle 
moving in a vacuum experiences a friction force proportional to its velocity 
as first discussed by Einstein [102,129,130,131,132]. Again, the friction force 
arises from the fluctuations of the electromagnetic field. A major difference be- 
tween friction forces and Casimir force is that the leading contributions come 
from the low frequencies for friction. By contrast, high frequencies contribu- 
tion gives the largest contribution to vacuum energy and therefore gives the 
leading contribution to the Casimir force. The friction force has been studied 
recently for a particle in close proximity to an interface. This force has been 
used to develop a near field imagery technique called shear-force microscopy 
[133,134,135]. The exact origin of these forces is still unclear. Recently, ex- 
periments on shear force [136,137] have been conducted in ultra-high vacuum 
in order to eliminate forces due to a water monolayer for instance. Differ- 
ent models have been developped to explain the origin of the friction force 
between two parallel surfaces [40,138,139,140,141] in the framework of fiuc- 
tuational electrodynamic fields. The case of a particle moving parallel to a 
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surface has been considered in [142,143,144]. The friction force for particles 
moving perpendicular to the surface has also been considered experimentally 
[145] and theoretically [146]. Volokitin and Persson [147,148] have shown that 
the friction force is much larger for the perpendicular case. The problem is still 
open because the discrepancy between theory and experiments is very large. 
Recently, it has been proposed [149] that the friction observed in the experi- 
ment by Stipe et al. could be due to the dielectric located below the gold film. 
Volokitin and Persson [147] have suggested that adsorbates might produce a 
very large enhancement of the friction forces. Clearly, more experiments are 
needed to clarify this issue. 



9 Concluding renicirks 



Many years after the discovery of surface polaritons, new discoveries and ef- 
fects are still being reported. A major reason is the development of near-field 
techniques that allows to probe the properties of surfaces with a nanometric 
resolution and motivates further work. Although thermal excitation of surface 
waves had been studied in the past, their major role in many phenomena has 
been realized only recently. It has been shown that heat transfer is dramat- 
ically enhanced in the near field due to the resonant contribution of surface 
waves. It has been shown that the electromagnetic field emitted by a thermal 
source may be partially temporally coherent and partially spatially coherent 
in the near field. It has also been shown that the far field emission properties of 
surfaces can be engineered by exciting and coupling efficiently surface waves. 
Partially spatially coherent thermal sources have been realized. Applications 
to emitting light devices arc under study. The role of surface excitations in the 
Casimir force has been clarified and the agreement between experiments and 
models is now satisfactory. The situation for friction forces is far less clear. 
There are still large discrepancies between pubhshed data and models. 
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A Green's tensor 



A.l Green's tensor in a vacuum 



The Green's tensor in a vacuum is defined by the relation: 
E(r) =//oa;'G(r,r')p. 



(A.l) 



Its exphcit form expression is given by: 
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(A.2) 



where R = |r — r'|, I is the identity tensor and u^u^ is a tensor in dydadic 
notation such that (u,. Uj.)A = Uj.(ur ■ A). 



A.2 Green's tensor above an interface 



For the plane interface system, it is convenient to use the representation due 
to Sipe [45] that consists in a decomposition over elementary plane waves. We 
use again the dyadic notation for the tensors. For instance, the s-component 
of the electric field is given by ssE = s(s • E). In the case of a Green tensor 
relating currents in the lower half space (r' in medium 2) to a field in the 
upper half space (r in medium 1), one has: 
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In the expression s — K x z/|i^r| and pf 
transmission factors are defined by: 
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Note that the transmission factor for p-polarization has a pole that corresponds 
to the surface wave. Thus the Green's tensor contains all the information on 
surface waves. 

In the case of two points lying above the interface in medium 1, the tensor 
can be cast in the form: 
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where the reflection factors are : 

p -ei72 + 6271 , 71 -72 

^2 = ; ^^12 = (A.8) 

ei72 + £271 7i + 72 



A. 3 Green's tensor for a two-interfaces system 



We now consider a two layers system. The upper medium is denoted 1 and lies 
above z — d. The lower medium {z < 0) is medium 2. Medium 3 is defined by 
d > z > 0. The Green's tensor relating the currents in medium 2 to the field 
in medium 1 with a film of medium 3 between media 1 and 2 are given by: 
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Note that the two-interfaces Green's tensor is very similar to the single in- 
terface one, except that the single interface transmission coefficient has to 
be replaced by a generalised transmission coefficient taking into account the 
multiple reflections. 
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When r and r' are in the film (medium 3), the Green's tensor can be cast in 
the form : 

^EE f (PK^EE 

G (r,r',^) = y^g (K,z,^')exp[zK(R-R')], (A.12) 
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where g (K, 2;, 2;') is the sum of three contributions : 
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B Fluctuation-dissipation theorem 



In this section, we derive the apropriate spectral density for an absorption 
measurement and for a quantum counter measurement. The operator needed 
to describe the absorption of a photon involves the normally ordered correla- 
tion function of the field operator < E^. ^Ei^^ > where E^^^~'^ is defined using 
only the positive or negative frequencies of the spectrum. We quote from [31] 
the cross-spectral density of the normally and antinormally ordered electric 
fields: 



£ii\r, r', cu) = ri{-uj)fXohu;^[l + coi\i{nw/2kBT)]Im[G^i^{r, r', cu)], (B.l) 
4f (r, r', uj) = r]{uj)iJo^^[l + coth{huj/2kBT)]Im[G^i^{r, r', a;)], (B.2) 

where rj{uj) is the Heaviside function. It is seen that the normally ordered cross- 
spectral density involves only negative frequencies whereas the antinormally 
ordered correlation involves only positive frequencies. They can be viewed as 
an analytic signal. Despite the apparent symmetry of the above equations, 
there is a critical difference as will be seen in the next subsection. 
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B.l Absorption measurement 



Let us now derive the time correlation function for the electric field. Prom [31], 
we have 



{Ek{r, t + T)Ei{r', t)) = 2Re[{Ei-\r, t + t)E\^^ (r', t) 



(B.3) 



Using the cross-spectral density (B.l), we get 



{Ek{v,t + T)Ei{v',t))^2Re 



[ ^ ex.p{-iu;T)tjLonu;'^Im[G^i] [1 + coth{fiu;/2kBT)] 



We note that 

/m[Gf,(r,r',-a;)] = -/m[Gf,(r, r', u;)]. 



(B.5) 



and 



[1 + coth(;k^/2A;BT)] = -2/[exp{h\uj\/kBT) - 1] 



(B.6) 



so that the correlation function can be cast in the form: 



{Ek{r,t + T)Ei{r\t))^2Re 



2^ 



ex.p{iu;T)2fjLou;Im[G^i^]e{u;, T) 



.(B.7) 



B.2 Quantum counter measurement 



We derive now the time correlation starting from the antinormally ordered 
correlation function. This choice is apropriate for a quantum counter exper- 
iment. The Casimir effect that depends on the total energy of the system 
pertains to this category. From [31], we have 



{Ek{v, t + T)Ei{r', t)) = 2i?e[(4+V, t + t)E\-' (r', t) 



(B. 



Using the cross-spectral density (B.l), we get 



{Ek{r,t + T)Ei{r',t)) = 2Re 



du 
2^ 
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We note that 



SO that the correlation function can be cast in the form 

duj 



{Ek{r,t + T)Ei{r',t)) ^2Re 



[ ^exp(-iu;T)2iiou;Im[G^i^][hu; + e(u;,T)] .(B.ll) 
Lo J 



We conclude by noting that both results have a similar structure and can be 
described by an effective spectrum defined for positive frequencies only. This 
analysis provides a justification for the heuristic argument often used to drop 
the vacuum energy fiuctuation. 
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